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TIC BENDING OF LAYERED PLATES — | 

Pister,! A. M. ASCE and S. B. Dong,? A. M. ASCE 


A system of equations governing the non- -linear bending of plates es 

“of two or more bonded, thin layers is presented. _ The effects of both trans- 


- verse and in-plane loading and thermal gradients through the plate thickness _ 


well as over the plate faces a1 are included. oo equations > 


al es enti in plate theory are employed. (See, for | example, 


Conventi 
Reference tf. symbols are defined below. 


oe oon | way temperature difference with respect to an initial 


ax. ne - linear coefficient of thermal nal expansion. 
- thermal : stress ‘and moment resultants 
- displacement vector ofa point on the top 


surface of the plate. 


xo» ? yor Vx strain tensor or components on | the top surface of of the 


- composite plate parameters 
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INTRODUCTION 


paper is with formulation on of a system of govern- 
ing the elastic bending of plates consisting of two or more thin, bonded el 
of isotropic materials. The equations are developed within the framework of 
< von Karman theory of bending, which ignores the effect of transverse 
shear deformation, but accounts for ¢ deflections which are of th the order of the s 
plate thickness. effect of transverse and in- -plane forces as well as 
_ thermal g gradients through the plate thickness and over the plate faces is a 
Although problems of layered plates have some attention on in 
struction. A general treatment of the linear theory of bimetallic shells and 
- axisymmetric bending of circular plates including thermal effects has been _ 
given by Grigolyuk. (2 j An independent formulation of the latter problem has 
_ also been given by Vinson, (3) who further restricted the problem to layers 
_ with the same Poisson’s ratio. Numerical results for a number of circular — 
= problems appear in a paper by Korolev. (4) The objective of the present 


Specification of the Temperature Distribution 


9 In many problems of interest coupling in the heat fi flow equations may be a 
"neglected and the temperature distribution in each layer determined from ap- 
- propriate solutions of the Fourier heat conduction equation. In —_— — 
in ea layer c can be expressed as 
Tx,y,2,t) = hy hy 
where a plate an with the associated snentionte system is shown in 
‘Fig. 1. Assuming no contact resistance at internal interfaces, the iaiaee- 
ture distribution will be continuous through the thickness. "Proper account Be 
contact 1 resistance can be taken with no difficulty. — For sufficiently thin a. 
the temperature distribution through the plate approaches piece-wise linearity. 
_ For time-dependent temperature distribution, the ne response | of of the plat 


= be ee as 8 quasi- -static in most cases. 
of the Large 
e 
the. coordinate system shown are: 
will 


where Yoo ' 


‘Karman theory are: 


— 
| 
Su 
| 
4 
nents of a po artial | 
Wo, d to be functions of x, iated with the von 


ing gives: 


ye 


vile ) “ten 


* 


2zw 


the n n plate laye layers. 
N. 
ox. ty)? 


xx 


-*Results obtained for beams indicate that this is a satisfactory approximation _ 
provided that temperature and/or — pressure Gradients in t the x, | 


directions are not severe. 5) 
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‘The thermal moment associated with the 


35 


If 


In the sequel the following additional constants will be 


% is convenient to regard the transverse displacement w and hin aan re- 
‘sultants Nx, Ny, , Nxy as primary dependent variables. . Thet three equations of | 
equilibrium for the plate element supplemented by the compatibility equation “_ 
_ for the in-plane strain components at the top surface of the plate constitute 

_the determinative | system for for these variables. These equations are: 
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are satisfied identically by introducing the Airy Stress . 


Substituting (20) in (17) )gives 
o* 74w- = — Mr] +9 + xy (21) 
and (22), from which w and F can be obtained, 
_ of the von Karman equations, modified to take into account the properties ee 


a plate made of bonded layers of different materials and to include the effect 
of temperature gradients.* The presence of the term KoV4 F in Eq. (21) an a 


= 


Ke V4 w in Eq. (22) indicates the coupling effect introduced by the layered con- 
struction, i.e., in-plane forces affect bending equilibrium. Since the four we 
boundary conditions associated with the system of equations do not differ 
those which obtain in the theory for single-layer plates, no further mention 
_ will be made here. In view of the difficulty in obtaining solutions of Eqs. (21) - 
and (22) even in the case of a homogeneous, single layer plate it seems desir- 


_ able to investigate the form taken by the in in particular cases. 


Reduction to Layer Homogeneous Plates 

For a layer Eqs. (21) and (22) take the 


ere D = v2)? N and MT represent arose and 
* 12(1- v2) 


7 -*Similar equations were obtained for a bimetallic shell of constant curvature a 
by Grigolyuk. (2) Unfortunately, aowever, the equations were not simplified to 
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E LASTIC BENDING 


modified to include thermal effects. the non-linear terms are 
neglected in Eq. (24), uncoupled, small deflection equations are obtained. Eq. 
(24) is then the thermo-elastic compatibility equation for plane stress and 
Eq. (23) is the equilibrium equation for a plate subjected to in-plane and 
transverse loads combined with thermal gradients. . Ifa linear temperature 
_ gradient through the plate thickness is assumed, these equations reduce to the 


_ standard form shown on pages 27 and 59 of Reference (6). ui ae a, 


for Constant Poisson’ Ratio 
oe In many instances the value of Poisson’ s ratio may be very nearly constant _ 
for materials with appreciably different moduli of elasticity. Assuming this 

GD be the case, the following semedaionaaia of the constants appearing in Eqs. 

( 


— 


Accordingly, these take the simpler form 
=v Ny - Mr] + +9 + xy — 


transposition of solutions and use of for n new problems. 


‘For deflections which are small compared to plate thickness, the quadratic 

7 terms in the strain-displacement Eqs. (3) are dropped. This results in set- 

7 ting equal to zero the corresponding terms in the compatibility Eq. (22); all - 
other equations remain the same. Accordingly, it is possible to simplify the as 7 
fundamental system as follows. From the above simplification, 

Eqs. (22") and (at) may now be combined to give the single Eq. 


i s in the small deflection theory the in-plane forces and transverse de- 


Equations of the Small ‘Theory 


flection are uncoupled, Eq. (21") for w can be solved in terms of prescribed 
in-plane forces and/or transverse loads and thermal gradients. | Note again 
that solutions may be obtained directly from corresponding single-layer plate 
4 probiems by transposition. The strain components at the top surface of the — 
¢ _ plate can be found in terms of the specified forces Nx, Nxy, Ny from Eqs. (18) =: 
the the deflection w obtained from the solution of Eq. (21" 
= of the one-dimensional large deflection problem of a homogeneous — 
strip subjected to temperature changes and transverse load been } 
given by Williams, 7 
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view of the invariance of the operators occurring in and (22) 
transformation to polar coordinates is limited to transforming the terms on “A 


s. (21) and (22) assume the pee 


ar 


where aj are constants. 
(8), (9) and (19) for stress and may | be 


transformed into a] appropriate expressions in polar coordinates. s. The elastic i 
constants will of course remain unchanged. f 
a Small deflection equations can be obtained from (26) and | (27) by dropping — 
aa -linear terms as for rectangular | plates. However, the problem is more + 
conveniently formulated in this case by choosing ug, the top surface radial 
displacement, and w, the transverse displacement, as primary variables. Sub- 
_ gtituting stress- resultant and stress-couple a, in the two —— 


equilibrium in polar coordinates leads to (8): 
[-DNy +CMy-CQ] 
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it will be nivel 1 om the equations are ‘not coupled and solutions for various — 
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Although transverse shear deformation is neglected in the ‘theory dis- 
" cussed, transverse shear resultants can be determined from equilibrium c con- — 


These stress resultants a: are determined from the moment equilibrium ‘Eqs.: 


7 


Combining these Eqs. with (19) gives. 


The distribution of shear stresses Tx2, T yz is e ob obtained by co con-— 
sidering first the interface shear stresses doncted by Tx’ and Ty(k) in wie 
F ig. 2. Equilibrium of the finite element in the x and y directions requires — 


=| 


the stress in the k th layer can be written 


S. “Theory of Plates and Shells”, McGraw- -Hill Book 


a Company, Inc., New York, 1940. 


. Grigolyuk, E. L, “Theory of Bimetallic Plates and Shells”, Akademiia = S 


Nauk, USSR, Inst. Mekhaniki, Inzhenernyi Sbornik, Vol. 17, pp. 69-120, x S&S 7 
Vinson, J. R., in Laminated Circular Plates”, Present al 
 edat the 3rd U. S. National Congress of Applied Mechanics, Brown Uni- a 


versity, Providence, Rhode Island, June 1958. 


4. Korolev, ¥..R “Thin Double- Pcie Plates and Shells”, Akademiia Nauk 
USSR, Inst. Mekhaniki, Sbornik, Vol. 22, Pp. 98- 110, 1966. 


sa) | 
ow @ 
1, 
| 
4 
21") 
asily 
ng 
(28) 
29) 
(28) 
ous 
4 
& 


Boley, B. A., “The Determination of Temperature, Stresses, 
* flections in Two-dimensional Thermoelastic Problems,” J. Aero. Sci., 
Vol. 23 (1), pp. 67-75, 1956. 


ion 
‘Subjected to Normal Pressure and Journal of Applied Mechanics, 

Dong, S. B., “Elastic ‘Bending of Layered Plates”, M. S. Thesis in ncivil 
Graduate Division, University of California, ‘Berkeley, 1958. 


4 


— 

— 

— 
_ 

| 

| 


October, 


ENGINEERING IANICS DIVISION 
Proceedings « of ‘the American Society of Civil Engineers 


SURFACES ON VISCOELASTIC /SUBGRADES@ 


B.C. floskin, and E. H. Lee2- 


‘ABSTRACT 
a ‘This report is concerned with the analysis of si stress and deformation of e.. 
subgrade strengthened by a flexible surface plate. subjected to a load distri- all if 
bution. Rate effects, creep and stress relaxation in components of such foun- : 
dations are included through the application of linear viscoelastic theory. Se 7 | 
shown that elastic analyses for such situations which have been 2 de- 
_ veloped in recent years can be directly extended to treat the analogous 
problem including viscoelastic effects. ‘Particular are presented 


“lutions, the rate influence causes variation of subgrade 


; and plate bending moments as the constant loads are maintained for extended — 
on periods. The significance | of the problem and the discussion of particular so- 
lutions are given in Sections 1 and 4 of the report. Sections 2 and 3 deal with © % 
_ the general development of the theory, 
[INTRODUCTION 
A problem of interest in in several fields of civil engineering 
that of determining the efficiency of a plate as a means of transferring a rela- 
tively localised applied load to a flexible foundation (Fig. 1). This «og 
arises, for instance, in the design of footings for columns and -_ ies meee 


Note: Discussion ope open until March a, 1960 0. Toe extend the closing month, 
__ written request must be filed with the Executive Secretary, ASCE. Paper 2195is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings — 
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7; a The results presented in this paper were obtained in the course of of See } 
Sater by the Office of Ordnance Research, Dew. of the Army, under 
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design of airport pavements. Most of the theoretical work aimed at obtaining 4 
rational design formulae for the above problem has dealt with elastic plates 
on elastic foundations and results of proven practical value have been obtained : 
on this basis. (1,2, 3) ) Nevertheless, it often happens that real foundations are 
‘not truly elastic; in  aeadinniee they may be subject to time effects as is gener-— 
ally the case for a soil foundation. | In this paper elastic plates supported by 
visco -elastic foundations are considered in order that analytical formulae, 
‘suitable for application to situations wherein time effects in a foundation aa 
_ of importance, can be obtained. Some preliminary work on this problem has 7 i 
already been reported(4) and recently the related of on visco- 
_In recent years it has been shown that in the stress. 
and deflection analysis of a system including linear visco-elastic components 
can be treated in terms of the analogous elastic problem having the same =i 
geometry and boundary conditions. In the case of a system, the geometry of | 4 
_ which does not change during the loading program, , and for which boundary _ 
conditions remain of the same type at each boundary point, for example pre- 
‘scribed load or prescribed deflection, application of the Laplace tri transform 


al 


= Loaded flexible surfaces over subgrades with surface 
ing fall into this c category. The present report shows how general problems of 

il this type can be treated by this method, and illustrates the procedure through — 

-, Linear viscoelastic materials comprise those for which the relationship a 
between stress and strain is a linear differential or integral operator with re- _ 
to time. . In addition to instantaneous elastic response to stress, such 
3 materials exhibit delayed elasticity and viscous flow, and for example 1 respond — 

a of constant as indicated in Fig. laws: are commonly 
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ars RESULTING 


STRAIN VARIATION 


& 
where o is the stress, the strata, isa the 
Pr Or , the being material constants and t the and Q 
of similar form: r Relation (1.1) is s written for ‘uniaxial ‘response - 
‘tw stress; in n the case of combined stresses analogous relations apply for 


stress components. For an isotropic viscoelastic material, two such laws are 
~ needed, one for shear effects and one for average hydrostatic tension and _— = B 


aS 


associated « dilatation. In general the response of such materials to stress can j 
be considered by replacing the elastic constants of an elastic body by the ee “a 
quotient of such rate operators. Relation (1.1) includes the behavior of ma- | 
terials represented by the / commonly used viscoelastic models of springs and > 
_ dashpots. Ih particular the Kelvin delayed elastic material and the Maxwell | 
material correspond to operators Q of first order, and operators P of zero i: 

and first order respectively. . The four element model, which is the simplest 

uf ‘relation to exhibit the three main 1 qualitative a: aspects of viscoelastic behavior— = 

4 instantaneous elastic response, delayed elastic response, and viscous flow as 

_ depicted in Fig. 2—is by quadratic operators. 
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4 
_ ‘Linear 1 visco- <ebnatie theory as expressed by (1. 1)o or by the more general 


_ of the } creep or relaxation function in the form of an n hereditary integral: 

: can be used to” reproduce | any y measured creep or relaxation law aw. However, — 

linear viscoelasticity implies the truth of the superposition principle, so that 


_ unless this is exhibited to a sufficient degree of accuracy by the properties of d 
the subgrade material under consideration, non-linear theory, and an inordi a a 


- nately more difficult analysis, , will be necessary. This constraint demands 7 
7 that if a creep curve  €o(t) is obtained at a certain stress So, , then the appli- # 
. - cation of twice the stress, 209, would have doubled the strain at each instant, 
that is, the new creep curve would be given by 2¢,(t). For many materials this + 
‘principle applies to a satisfactory degree of approximation as long as the ae: z 
stress does not exceed a certain critical value. (7) presents a review of stress _ 
analysis in linear viscoelastic and of methods of measuring materi-— 
ease of application of the method mentioned above, of transforming 
_ problems of visco-elastic stress analysis into analogous elastic problems de- 7 q 
_ pends on the manner in which the elastic constants appear in the final ex-_ t 
_ pressions for stress components. If this is in the form of simple rational — 1 S, 
functions, then the evaluation of the equivalent visco-elastic problem is ; ae 
straightforward, , but if the elastic constants appear as arguments in 
cendental functions, difficulties of interpretation arise. _ For some problems rl | 
there may be a choice of representation of the solution which will provide — 
appreciable simplification from this standpoint. In the next section some = 
elastic foundation problems are discussed in a form which is suitable for ex- Fi 


We will be concerned with axially symmetric loading in 
" dations for which edge effects can be disregarded, so that the foundation and - 
14 surface are effectively i infinite in extent. We will consider loading problems > 
_ in which inertia forces due to the deformation are negligible so that quasi- pin 
. static analysis will be applied. This will introduce no restriction for most mM < 

foundation problems, although there will be a limiting speed for a moving 
vehicle, or for the build-up of a dynamically applied load beyond which inertia 
‘forces will have an appreciable influence on subgrade stresses. = ~ ale 
_ For axially symmetrical problems both in three dimensional elasticity j 

_ theory, and in plate theory, stress components which vary radially in ‘pro- a 
portion with the Bessel function J)(mr), where r is the radiusandma 
oe 4 constant, play a particularly significant role. For elasticity theory in cylindri- 

i 4 cal coordinates, this distribution produces separation of the axial coordinate _ 


z, and the radial coordinate r in the basic equations and so generates simple 3 
solutions i in closed form. For example, for a load distribution q(r) = =Jo(mr) 
on the surface of a semi- half-space, the defiection 
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4 al layers of different « ro media ¢ can be treated in a pos Mint lone and that 
_ the surface deflection is given by an expression analogous to (2.1), with the © - 
coefficient of the —— function a more complicated, but rational function = 4 
— The analysis for a single dinette layer resting on an elastic alf- -space of 
gaa material properties is more simply treated if the dimensions of the 


permit thin plate theory to be used to the layer. 


where D is the plate flexural rigidity, s(r) is the reactive pressure senare saute’ to 
~ the plate’s lower surface by the foundation, and q(r) is the intensity of the 2 
plied load. Owing to the radial symmetry present in the class of problems 


under consideration, ,W, 's and ¢ q are functions only of the radial coordinate r, i 
and the operator is given by 


- om 23 (mr) 
of (2.1) with (2. and the identity (2.3), gives the deflection 
i; 7 the plate on the elastic half space for the surface loading q(r) = Jo(mr): | ‘ : 


A similar result is obtained if the elastic half space is sili by the the S 
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_where k is a material constant for the foundation. Such a a ae is ‘a 
commonly termed the dense liquid type, since it represents the flotation eat 2 

t 2 a dense liquid would i produce, although the value of k ‘ normally associated with a 
_ foundation problems: is much larger than the e density of the material displaced. = 
Most of the formulae used for the rational design of conorste pavements have th ; 
been obtained by assuming this type of foundation response. For this case Me 
(2. 2), 3), and (2.5) determine the deflection of the the load 


Reissner(12) ) has d a theory of subgrade action which i is ir inter 


mediate between the above two assumptions. It is — the 
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FIG. SMPLE ELASTIC FOUNDATION 
_ theory, and supplies some of the ‘shear coupling which it is 3 absent in the liquid t 
‘subgrade assumption. 
- contain space derivatives only in the form of V2 in the surface, the use of the he 
Bessel function pressure distribution \ will- achieve Separation of the ‘space and 
_ Thus for all the foundation representations considered, the application of | 
Bessel function loading distribution leads to simple deformation ex- 
pressions, and incidentally, expressions in which the elastic constants appear ; 
as simple, multiplying factors; these lend themselves readily to generalization 
for visco- ~elastic ‘materials. However this load function introduces certain — he 
difficulties. eH (mr) represents a a damped oscillatory load with maximum in- q 
tensity at the origin. The amplitude of the oscillations decrease to zero with 


increasing radius, and in fact the function is represented 


% 
This ensures that the plate deflection and its derivatives approach z zero as sj 
increases, so that all moments and stress components do likewise. However, = 
the total load on the surface is not he total load to b 


8) 


As b increases this oscillates with an amplitude that t approaches wen so 


— example, a distribution of load that is of uniform intensity Pagpcade an area 
of radius a, but which is zero outside this area (as depicted in Fig. 1) canbe | 
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EM4 SUBGRADE SURFACES 
| Since r appears only in the e term Je (mr), the analysis detailed above for de-— 
flection of a flexible surface ‘supported by certain types of foundation, can be — 
applied directly to the integrand of (2.9), resulting in an infinite integral over e 
_m for the deflection and other quantities which may be of interest. For ex-_ 


ar ple, the deflection ofa plate resting on an half- for load- 


a 
am 


ing is given from (2.4) by the 
of this type is carried out in the 
In spite of the difficulty of not being able to define the total load on an » infi- 
nits plate for the load distribution Jp (mr), it may still be useful to use a finite 
sum of such terms to represent prescribed radially symmetric load distri- Bae 
_ butions adjacent to the origin. Again the total load for the infinite surface will 
“not be defined, but since the deformation under the load is effectively a local — /? 


phenomenon, and the stresses at infinity approach zero, it is reasonable to ex- 

pect that if the series of Bessel terms represents the load distribution satis- 

-factorily around the origin, it will be a usable representation for the determi-_ 

nation of deflections and stresses in this neighborhood. The reason — 

suggesting this, when an integral representation of the type (2.9) is —— a 

which avoids these complications, is associated with the simplicity of evaluat- 

ing a finite sum rather than an infinite integral. If an elastic layer analysis - 

| type developed by Burmister(3) is ‘used, the complexity of the displacement 

| expression for a single di (mr) loading term n may dictate the use of a finite ; 

Series to achieve a solution within practicable computing limits. 

_ An expression of the type (2.9) is developed from the orthogonality of 
Bessel functions in a manner ‘similar to the development of the Fourier inte- Ps 
gral for trigonometric functions. Corresponding to the Fourier series, watch te 
represents an arbitrary function over a finite interval, the Fourier- -Bessel 

series of Bessel functions also achieves this. It is developed in an analogous 
manner (see for example (9) Chap. VIII) by means of orthogonality relations a 
x over a finite interval. _ Whereas the Fourier series represents a periodic | aol 

. function which ‘repeats itself outside the chosen interval, the Fourier-Bessel 

_ series in terms of Jo(mr) also introduces a continuing contribution with in- _ 

- creasing radius, but this is attenuated as the radius increases and approaches — 
zero for large radius as does the Bessel function itself. ‘Fig. 4 ‘shows the 

; _ representation by a finite number of terms of a Fourier- ‘Bessel expansion of 

;, a uniform distribution of pressure of the type depicted in Fig. lof magnitude ant 
2.5 over a circle of radius 0.2. Fig. 4a shows the result of summing five ae 
terms, and Fig. 4b, thirty terms. In each figure the broken line shows the 7 : 
distribution to which the finite series isa an 1 approximation. The! 


= figures show that these series introduce a ring of negative pressure <a 
 houanenriens of r = 2, and such rings of decreasing amplitude will ap- 
pear for larger values of the radius. _ These distributions suggest that the brs: 
finite series will provide a satisfactory representation of the single region of 
‘alors pressure, as long as the deflection solutions associated with the a 
rate rings of pressure introduced by the finite series do not interact appreci- 
ably. This can readily be checked by computing the deflection distributions. — 
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. prec for a ‘limited time duration, after which the deflection distribution 
- spreads out and introduces interaction between the separate pressure regions. 
The finite series representations depicted in Fig. 4 were based on the interval 

r<l. A interval wouid displace the first negative pressure ring 
2 to a larger radius, ‘but would call for: a larger number of terms to achieve — "on 


the accuracy ‘needed for the approximation to the desired load distribution, it 
should be born in mind, as shown by Reissner, (2) that in the range of loading 
areas and surface plate thicknesses for which the subgraded pressure is ap- ¥ 
ft preciably less than the load intensity, that is the practical design range, , the q 
deflection and subgrade pressure functions are insensitive to the detail of the 
form of the load distribution. Thus the approximation shown in Fig. 4ais 
likely to be quite adequate for. ‘computing foundation properties. Although for my 
_ the particular problem solved in detail in this report numerical integration _ 
™ using a Fourier-Bessel integral relation of the type (2.9) was utilized, aspects 
4 of the application of finite series have been discussed since, as mentioned _ 7 


~ above, the simplification achieved may b be worthwhile when more complicated 


3. Analysis for Viscoelastic Materials aft 
a asi ‘The previous section dealt with stress and deflection solutions for elastic 
flexible surfaces and subgrades. _ When there are components which respond > 
to stress in a viscoelastic manner, , corresponding laws must be introduced, 
=: the resulting stresses and deflections will be found to vary with time even 
_ When the applied loads remain constant. The shape of the distributions Te 
stress and deformation will in general also vary. A convenient method of a 
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where | oh isa neomee: el it is well ale (see for example (10) Chap. 1) = 


the Laplace transform of the time derivative of a function is simply p times 
the on transform of as function itself if the function is zero initially: 


with this viscoelastic | laws of the form a. become 
i equivalent to elastic laws with the elastic constants in the form of rational 
: functions of the parameter p. Thus the determination of transformed varia- i 
bles, of say stress or deflection, becomes an elastic problem of the type dis- 
= cussed in the previous section, and standard methods of inversion of the Tae 


' _ Laplace transform | determine the stresses and deflections as functions of the — 


onsider, for example, Eq. (2.1) for the deflection of the surface of an wd 
elastic half-space subjected to the load distribution q(r) = Jgo(mr). Let ns 
‘suppose the material to be viscoelastic with the shear modulus G reptnced wa 


the viscoelastic 
where T and y are shear stress and shear strain respectively. Further, let 
us suppose that the operator for the visco-elastic body which corresponds to er. 
-Poisson’s ratio » for the elastic material, is a constant, which means that in a 
a simple tension test of the viscoelastic material, the ratio of lateral to — . 
tudinal strain remains constant throughout the test. _ This would be the case, 
for example, were the material incompressible, corresponding to » = 1/2. a. a. 
Fora viscoelastic problem, load loads varying with time met t usually be be 


, and so let us take: th 


ion 
ns. 
val 
i. 
| 
— 


2 
It is convenient for analysis, and is called proportional loading, since at any 
two given times, tie load intensity at all points changes in the same ratio. 
with the Laplace | transform on (3.3) and (3.4) 


= 


In terms of the transformed variables, the the deflection vie the ani os 
viscoelastic solid corresponding to (2.1) 


by replacing in (2. 1) elastic. quantities ie their’ viscoelastic equivalents. s. This 
step is based on the fact that for transformed variables the viscoelastic oa ie 
problem | becomes an elastic one. in the form: 


9¢ 
and bearing in mind the linearity of all 
functions of space variables can be taken outside the rate operators, we see 
that there is a direct correspondence between the shear strain response fon 
the applied stress variation for uniform shear, and the deflection response 
w(r, t) to the time varying load $(t t) Jo , (mr) for any \ value of the radius of ee 
surface loading problem. One could, for example, ‘measure bee shear strain ‘ 
fora certain history of stress application in uniform shear, and deduce the 
‘surface deflection for a proportional variation of surface loading. In this ‘ce 
the additional factor which arises is expressed i in terms of » and m, which 
will be known constants for a specific problem. In this application only the 
‘single operator | quotient Q/P appears as a multiplying factor, so that the 
FS lution can be determined from direct measurement of a related situation under 
_ homogeneous stress conditions. hs This is even true for arbitrary radially sym- 
‘metric proportional loading on the surface of a 1 homogeneous half space. Such | 
an arbitrary can be represented by the Fourier-Bessel 


‘on fu unction A(m) determines the shape of the load | distribution, 
be o(t) the: variation with time which must have the same form at all surface — 
Sav _ Eq. (2.9) gave an integral of this type for a cylindrical distribution coll 
uniform pressure of radius a and intensity qo, for which A(m) = = qoaJj(ma), — 
and no variation with time was considered. Using the Laplace 1 transform pro- 
_ cedure, and the solution of the associated elastic problem in the form (2.1), i 


obtain for t = ofa half- — the of 
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elastic G. The term - ty u ( ) can be the 


because it does not contain m, the integration variable. Applying the inverse a 


yt) = vt) 5 4 


The combined operator (1+, 6 t))/G(8/6t) cannot now be directly related 


— toa simple loading problem in homogeneous stress as was the case for 


_ constant .: and in general it will have to be treated for specific operators — aj 


4 


G= = Q/P and yu. ae If these correspond to the common models of springs and 
dashpots, which implies that the operators are rational functions of (d/at), 
_ the method of partial fractions can be used to deduce w(r, t) from ¢(t). . Simple 
examples of this procedure are given in the next section. 
-" It will be observed that in (3.9) the rate operators can be combined into = 
single operator which is a factor of one side of the equation. This means that | 
5) | to each particular shape of load distribution determined by A(m), a unique 
shape of surface displacement will correspond, being proportional 
ip? a shape will be retained throughout the motion, so that points at different — 


radii will be displaced proportionately, and this will take place independently 
of the particular form of the viscoelastic operators or the force variation 4 a 
| g(t). This situation is extremely special for visco-elastic deformation _ 3 


av _ problems for which in general the space and time distribution of stress and 

um _ strain both differ in form from each other and at different points in the body. — a 
pal, For example, a constant applied load can produce differing stress =. 


If we now consider a viscoelastic foundation wits a flexible surface, we a 

consider the transformed equivalent of (2.4): th ~ 
= (mr) (3. 11) 


{ 
te toad distribution is given 


8) i The plate flexural rigidity D is considered as a function of p to permit visco- 
id - surface layers to be treated, as well as viscoelastic subgrades. This 
ee may be needed to provide a satisfactory analysis for chemically treated as ae: 


Le 4 surface soil layers. In this case the rate operators cannot be separated out * 
of as a factor apart from the parameter m, even for constant D, so that when an =e 

arbitrary load distribution _ e is considered, the transform of the resulting RE 


* A 
| 
i 


- ‘The combination of the operators and m in the integrand indicates that for a 
fixed shape of load distribution, the shape of the displacement at 
will in general vary with time. A similar situation arises if 


__‘The combination of the application of the Laplace traneform ant the so- 
’ = of associated elastic problems already used in foundation analysis and f~ 
: _ reviewed briefly in the previous section, permits the stress and deformation a 
analysis of a variety of foundation problems involving viscoelastic components. f 
The complexity of the analysis depends on the particular operators needed wie 
_ represent the materials, and in general this increases markedly with the 


order of the if Burmister’s(3) elastic solutions are used for 


the previous section. gy ‘examples of the method are detailed below. 
_ It will be noted that the interpretation of elastic solutions to generate vieco- 
‘elastic solutions depends on the elastic constants ‘appearing ina simple nel 


a for all cases cited above, for which the elastic constants appear : asa 1 factor of 


the Bessel function with radius in the argument. For the evaluation of elastic — 


7 solutions, manipulation has been successful in avoiding integral expressions, - 


for example, the deflection of a plate on a heavy liquid type subgrade with a 
cylindrical region of uniform pressure of radius a and intensity qo can be ex- 
Pressed in closed form using the Hankel functions as follows: (2) 


m over the range 0 to ~ with the spectral function (ma). Clearly (3. 14) 

_ would be difficult to interpret as a viscoelastic problem with the elastic ee. 7 

constants appearing in such a complicated fashion in the Bessel function argu-— 

‘ments. Thus in developing elastic solutions for viscoelastic applications _ saves a 
‘special considerations arise which may dictate the suitable form of form of solution — 

different from that for most convenient elastic analysis. 


consider a subgrade covered by a flexible elastic plate, initially undis- 
turbed. At time t t = 0, a load distribution is suddenly applied, and 
it remains constant. Thus > 


purely 


by the Maxwell viscoelastic law. 


(Fig. 
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Spatial whe: 
bgrade, with its response to load 
= _ Thus each of the springs shown in 


the notation of Eq. (2.2), A and B being subgrade 
— law represents a material which exhibits instantaneous elasticity — 


combined with viscous flow if the load is maintained as is apparent from the 
_ model representation, Fig. 5. Under constant — * ‘such a material 
: _ The Laplace transform of the plate deformation Eq. (2. oT for the load 5 
_ distribution ( (4.1) and zero initial deflection is: ri wae 


transfor of (4. “4 


nd combining these we have the associated 


We consider the spatial ie al of applied load to be of uniform tmtenatty 
Io over a circular region of radius a about the origin and zero elsewhere. ae 
Representing this by the integral of Bessel functions (2.9), the solution of 
(4, aia is given by an integral over solutions of the type (2. 6): Ah, awe 
= q, 2 dm 
is expression can be inverted in an elemen ary fashion using partial — 
This: be ted lementary fashic arti na 


fractions and on performing this it is found that 
corresponding fot reactive pre pressure can an be most conv convenien 
iy by substituting (4.5) in (4.2) and then inverting. The result is | 


standard formulae of thin plate theory, the time being merely a i for 
this determination. . For bending moment unit 


m(r) (r, wie is given by 
sae 
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FIG.6. DEFLECTION PROFILE FOR PLATE ON MAXWELL 


ae where use ei aa been made e of the B Bessel function 1 relations patie maton 


_ These integrals were evaluated numerically by means of Simpsons rule. . The 
mé in the denominator of (4.8) calls for a limit evaluation for small m. For 

large m, the integrals converge r rapidly because of the properties of 

_ functions, and the inverse powers of m appearing. If the convergence for -_ ; 

: large m ist not sufficiently rapid, the well known ‘asymptotic expansions for tee 
_ Bessel functions can be employed to provide tt ain which can be evaluated — 

closed form a particular value value of m. 


/ 

Te 
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Profiles of the deflection and foundation reactive er: at various times 


3.7. REACTIVE PRESSURES FOR PLATE ON MAXWELL FOUNDATION 


region, i.e. p = 0, is shown in Fig. 8 for this same system. The particular 

was chosen for the parameter since, as pointed out by 

 Reissner, 2) this represents the approximate limit of effectiveness of the 
‘ flexible surface. As shown in Fig. 7, the maximum subgrade pressure induced 
is about 30% of the applied load intensity. For larger | values of A, that is ed 
re- 


either a larger area of applied load or a thinner surface plate, worthwhile 
duction of the pressure on the subgrade is not achieved. The subgrade ae 
= _—— decreases with decreasing A for a fixed load intensity, but ——— 
A < <1, the analysis becomes unreliable because the region of loading is too 
_ small compared with the plate thickness to permit plate theory to be applica- 7 
For small A, Burmister’s solution(3) based on three dimensional elastic 
would have to be employed. 
_ The influence of the viscoelastic property of the subgrade is clearly evident 
in Pies 6, 7 and 8 in comparison with elastic > subgrade analysis, which, for if 


av the constant load distribution applied, would predict a constant deflection and sf 


constant subgrade pressure distribution. Fig. 6 shows that the deflection 


BS pattern gradually spreads out, for while the central deflection doubles in the an 
interval from T = 0 to 4, the deflection at p=4 increases manyfold. Under 
the more intense central peak of subgrade pressure, the deflection at the 
centre of the plate increases rapidly, causing an increase in plate curvature ai 
«there, which is reflected in the increasing central bending 
‘moment shown in 8. 
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FIG. 8. BEN AT CENTER OF PLATE ‘ON MAXWELL 
7 

: the dashpots have had time to deform. This determines the subgrade pressure 
TP . ‘distribution denoted by 7 = 0 in Fig. 7. Since 7 represents dimensionless = 
mes | time in terms of the relaxation time, to = A/B, of the Maxwell type subgrade, — 
fort << tg the situation is closely approximated by the elastic solution. Thus” 


ed 4 for large i elastic analysis will be satisfactory for moderate loading times. 
2 As Tr increases, the increased flow of the central part of the subgrade relieves 

ee : the pressure there and the subgrade pressure distribution spreads out. The _ . 
_— applied load is thus supported over an increasing area of subgrade at reduced q 
uced Pay _ The effective increase in span of the plate carrying the load causes - 7 
the increased bending moment mentioned 
With increasing time, this effect will continue and will approach the analo-_ a a 
aaa gous situation to the floating beam mentioned by Freudenthal and Lorsch,(5) — 7 


ea _ When in the limit the subgrade pressure will approach zero, as a large area of 


00 _ the plate transmits deformation to the foundation at a uniform rate which ap- sx 
ca- | proaches zero. In practice this limit will not be reached since no practical — == 4 4 
stic _ subgrade system can permit unrestricted flow as does the Maxwell body, and ares j &§ 
«moreover, plate edge effects may become important as the deformed area 


ident 7 = _ In order to represent a practical foundation, the material properties af 
oe 9 mast be changed so that the deformation is in effect restricted by a spring in in = 
an with the dashpot to transform the viscous flow to delayed elastic ey o4 
a deformation. However, in such a case, the solution detailed above will _ aoe 
ae ; vide a good approximation for a limited time of load application. is 
Ay 
| Consider now material properties which correspond to restricting the flow 
+e of the Maxwell type foundation by, in effect, inserting a spring across the a 


re dashpot and slamenta of Fig. ure the viscous fic flow to elasticity. The 


1950 


; linear solid, and Freudenthal and Lorsch(5) gave an example of how this can 

approximate soil behaviour. The stress deflection relation for this foundation 


where A Bj and Cj are subgrade constants. Using this relation in place of b 
(4. 4.2) of the previous problem leads to a solution by an iw 
and the the final al expression for plate deflection corresponding to (4.5) is 


+1) (By Dm 


— 
‘rep 
sin 
syn 
oft 
aD 
det 
4 
4 q 
— 


it has been shown in this report that a variety of elastic stress and defor- — 
oo distribution analyses which have been used in the consideration of oy 
foundation problems can be extended to materials with a time dependent re- 
“sponse to stress represented by linear viscoelastic behaviour. ‘Thus if rate a 
effects, creep and relaxation behaviour of foundation and subgrade materials — ) 
can be satisfactorily represented by linear viscoelastic laws, analyses are 

- available to include these effects in the consideration of foundation design — 
problems. It may be that if linear viscoelasticity does not represent these = 
fects in foundations and subgrades quantitatively, qualitative information de- ; 


‘rived from such analyses will be valuable to supplement information derived 
a elastic solutions, from which such rate influences are entirely absent. 
= _ Throughout the analyses it has } been assumed that ‘no separation ee | 

- the ‘flexible surface and the ‘subgrade occurs, even though negative subgrade 
| pressures are exhibited by certain solutions. In general this is unlikely to be ~ 
|: significant difficulty, since the complete solution involves the superposition 


an of a uniform pressure due to the weight of the foundation m materials, and this - 

> will tend to eliminate the negative reaction pressures. Because of the lineari- 
ty of the analysis, such superposition is allowable without violating any of the 
equations used. Another aspect of this ¢ question which has been discussed by 

Holl, (11) is that the oscillatory solution with negative subgrade pressures oc- a 
curs with the heavy liquid type subgrade behaviour which is less likely to +a 
represent detailed aspects of foundation action than the homogeneous material 


subgrade. The latter does not lead to oscillatory ‘variations of displacement — 


Reference (1) gives s solutions for finite plates on subgrades u using 


“similar Fourier-Bessel series to those treated in this report. Also non- 7 
symmetrical loading is considered by expanding in series of Bessel functions 
of higher order of the radius and trigonometric functions of the angle in polar — a 
coordinates. _ The results are much more elaborate than those presented here, 
but in principle the present viscoelastic method of analysis can be applied ‘~ = 
without change. Individual cases would need to be studied to assess whether 
‘the computing task involved makes the final numerical evaluation feasible. ety : 
i. The particular « examples treated in this report analyse loads suddenly ap- i. 
plied at t = 0, and held constant thereafter. Varying loads can be treated in a 
ton manner by modifying (4.1) to include a function of t other than the step 


= 


function, and by including the corresponding transform in . (4-3) > Alternatively 
the solution already obtained for step function loading can be used directly in 
a Duhamel integral to give results for varying load. For example, any quantity 
"determined for the steady load (4.1), can be generalized for an applied load 
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> This relation arises from the wien of the analysis, and applies for othe 


quantities such as stress and bending moment. ar 
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Proc. SOLUTIONS FOR ONE-DIMENSIONAL STRUCTURAL LATTICES 
wit Donald L. Dean,! M. M. ASCE and Selmo Tauber2 

A brief t is presented on the a the advantages of using ‘closed form so- 
eam ee to difference equations as a tool for structural analysis. In solving 
bee "directly the difference equations obtained by substituting standard force- -joint 
he. _ deformation relations into the general joint equilibrium equations, ‘complete = 
4, | solutions are obtained for for regular parallel chord Vierendeel and triangulated — . 
jgesum ad it seems that methods familiar to persons conversant with the calculus of Shae 
a @ finite  differences(1) are seldom used as an analytical tool by stress analysts =e 
and structural engineers practicing today. Although in recent years this 
wwe, branch of mathematics has been applied most extensively by statisticians, it . ; 
eed ly is interesting to note that one of the near- -classics(2) on the subject deals pri- _- 
marily with applications to the field of structural analysis. Portions of 
r. for a work have appeared in recent texts,(3,4) but there have been very few a oe 
i published attempts to extend it. This is unfortunate, as modern operational _ 

Pa, methods in the calculus of finite differences may readily be used to obtain a 
rare complete practical solutions to a rather broad class of interesting structural _ > 
al _ problems. This is especially true in the area of structural lattices, i.e. wots fa 7 

in regard to loading can be dealt with An 
: — load distribution is selected for the Problems | here with the hope that they will 7 
Note: Discussion open n until March 1, 1960. To extend the closing date one month, a 
os _§ written request must be filed with the Executive Secretary, ASCE. Paper 2196 is mv a 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 


of the American of Civil Engineers, Vol. 85, No. EM 4, October, 1959. gon 


Associate Prof., , Dept. of Civ. Eng., Wale. of Kansas, Lawrence, Kans. abot a 
2 Assistant Prof., Dept. of Mathematics, Univ. of Kansas, Lawrence, cane 
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october, 1089, 


- thereby serve as henstettn samples of the approach. Load variations exp expressi- 
Y.. ble as algebraic polynomials can usually be solved and, at least in the case of © 
systems represented by difference equations with constant coefficients, tran-— 4 he 
gcendental variations in the loading are tractable. The step function (loading 
" _ begins at an arbitrary location) and the delta distribution (a single load arbi- — “a 
_ trarily located) are examples of the latter. The dynamic analysis of structur- 
- al lattices also can be accomplished by direct solution of the governing differ- 
~ The usefulness of their direct solutions in the design and/or ‘analysis of ; 
lattice structures leads the writers to recommend that consideration be given ot 
m to the inclusion of a study of difference equations in the undergraduate cur- 
--riculum of students | ‘specializing in structural design. their knowledge, 
such a program does not exist currently - The modern curriculum regularly | 
includes a study of differential equations, which area limiting case of ‘differ-_ ; 
ence equations. The methods for solving the more general case involve ida 
2 enough analogies to the approach used in differential equations so that one or 7 
= additional lectures will enable the student to solve elementary oneal i 
~ equations. Preliminary experiments with the introduction of such material 
a second course in indeterminate structures have been successful. _ Since 
in the first course he has learned a number of methods convenient for solving 
7 _ the coarse lattice (less than 20 joints or units), the student readily grasps the _ 
_ significance of a tool with which he can handle finer structural lattices. The i 
simpler solutions for regular lattices are often 


& 


In the two examples given here, a general solution is obtained which con- 7 tw 
tains enough arbitrary constants to satisfy any given set of boundary con- oe : 
ditions. For the Vierendeel truss, the constants number three for 6x and four 
for , aa x: This is a generalization of the results presented in reference (2) 
a where the shear rather than the load was used as a starting point. The so- ms 
- lution for example two seems to be new. Only energy due to bending is con- i, ; 
sidered for the Vierendeel truss and direct or primary stresses are the sole ~ 

consideration in the case of the truss. The loading is limited 


statically applied joint forces. 


_ An unknown deformation (slope-deflection) method i is used. The adi- mtidng 
mensional unknowns are: @,x, rotation of joint x and deflection of 
in x divided by h (see Fig. 1 at end of paper). The effects of axial forces os. 
on curvature are ignored so that typical moments in te terms of unknown defor-— 


The necessary di! difference equations : are obtained by substituting these re- 
ue a lations into the two equilibrium equations expressing: : 1. The summation of 
moments on x zero of interal forces on 
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Chord Sti = El/h)chord; 


Ratio, = Ky/Ke; 
es atio, = 
Intensity Constant, Co = 
Adimensional Load Variation Function 


Eqs. (1) and (2) are transformed to a standard terme by introducing the — 7 a 


- Boole displacement operator E, where E f(x) = = — and E™ f(x) = f(x+m), 

‘The we now be made more re specific by selecting Px=lie.a 

concentrated lateral load of p on each joint. _ (1a) and (2a) are then reduced to > 


two etciemaace equations by using determinants, as with algebraic equations. 


+1 
om 


in the following for reference. The given difference 


¥(E) isa polynominal of E, has a solution, 


| where r are simple roots of and “a” is not a a root of The 
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e of 33 
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zee a” isa 2 double root of ¥. if rj is equal to r2, the first saints second centr | 


_ Where B and 1/B are the roots of ¢. -, ’ 


‘The in n and (20), Cc, Cg, - - C4, represent only four 


independent constants as may be determined by substituting the ‘solutions into 
rors and sad J. In matching coefficients of Ii am terms, | it is see 


many instances, depending upon the type of boundary to 


-.evaluate the constants, it is more convenient to express the solution in terms 
_ of hyperbolic rather than exponential functions. Proceeding in this smal 


problem considered, one obtains 


a 


‘4. 


— 
| 
= 
= pax (1d) 
ng 
— 
a (1d), (2d), (3), and (4) comprise ace The constants Cy, Co, 
—— ierendeel truss un ituting the given boundary conditions in Par 


will n be evaluated for a specific set of boundary con- 

ee ditions to complete the example and furnish working formulas. The mixed 


ay boundary conditions of a cantilevered structure will be used (see (see ‘Fig. 1 at end — 7 q 


These conditions substituted into (1d)... (4) 3 expressions for Cy 
‘this specific case then, Eqs. (1d) . (4) can be written: i 
4 


4] 1- 3ex(n=-x) C {(en+1)tanh = sinh An or 


=| 


-2x -(2n+1)(cosh Ax -1) + sinh A 


ay 


4 


- (2d) A convenient check may be made by comparing results with those obtained 
=a by one of the “coarse lattice” methods for small values of n. tie Raat 
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equations expressions give the ‘member — in terms of 
_ the unknown joint deformations. Although readily available, (5) these cnnatiene 
are not so well known as the slope deflection equations and will be explained — 
_ in more detail. For either plane or space trusses, depending upon whether _ 
_ the indices assume two or three values, the equations wr written in im-compacted = 


Where is tote Fig. 2 at end of paper): the indices i and j indicate components _ 
parallel to one of the axes of the reference frame, here assuming values 1.and 
a (reader, Sa accustomed to the tensor notation, where a repeated ii index denotes 
summation, will find the first summation redundant). k indicates the 


7 


for joint k. (For convenience in writing the difference equations ae a the 7 
subscript k is replaced by a subscript and a variable, for example, the hori- 
zontal movement of the grd joint on the top chord is written ut (3) and the yt 
vertical deformation of the rth ‘N joint on the bottom chord is s denoted Up(r).) a 
7 The member { forces and and force Components, in terms of the joint de it defor- 


Here, the four difference equations needed to solve for ‘unknowns SUR), 


U, (r), ( tr) and (x) are obtained by substituting the indicated physical 


quantities into Eqs. for the : r™ top and bottom joints. 
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; p-1(r) and p~2(r) are adimensional functions of (r) describing variations Ca 
eee and vertical components of top chord joint loads and q-l(r) ar and Asttit 


ci. 2(r) are similar functions for the lower | chord joint loads. For example, 


The sample | problem 1 here will be the regular triangulated truss with uni- 
for form vertical joint loads at the upper chord level. Therefore, 
“Then Bas. (0), (8), (9), ont (10), operator notation, become: 2 


B-1)ul(r) + 0 + 
+ a) (r) + (r) 2u opm 
ie As in example one, these equations are subjected to the usual algebraic _ ‘ 
- operations — were used here) to yield four equations each with | 
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Pe _ The solutions a (7c), (8c), (9c), and (10c) may be written as algebraic poly- 
‘nominals of the third and fourth n degrees c containing three o: or four arbitrary 
constants ; respectively y depending « on the order of the gc governing difference ied 
equation. In this particular case, however, it is more convenient to work with > 
factorials and the final result is just as easily used numerically as if ordinary 

~ polynominals had been used. A factorial of the mth degree is written (rm 

and is equal to r(r-1)(r-2) . (r-m+2)(r-m+1). If. it is desired to transform 
a result which appears as a polynominal of factorials to its equivalent as = 

algebraic polynominal, a table of coefficients called Stirling (1) numbers is | 

used. The convenience of factorials is understood if one remembers that — “| i 
difference equations are sometimes more easily used if written with epleuiaes: 
other than E. For example, Eqs. (7c)-(10c) are more compactly written ae 
the forward difference operator A where Af(x) = f(x+1) - f(x) or A=(E-1). 
An advantage | of transforming difference equations from the E to the ‘A notation 
is that one may then use the difference formula = (m)p(r)m-n which is 
analogous to the well known derivative 1 formula for an operation on the power — 
form in infinitesimal calculus DA(r)m = (m)p(r)m-n The to 


5 g(t) + + + tag) (rel) 5-627 


which apparently could have been grouped with similar terms in the comple- sf 

_mentary solution by a simple change of constants. . This was done for algebraic 

convenience in obtaining relations between the arbitrary constants. _ Also, one — 

_ must exercise care in selecting particular solutions to the independent — 
equations, (Tc)- (10c), to avoid the difficulty of terms which will not satisfy ‘yy, 
the dependent Eq. (7a)-(10a), 
Inspection of the four independent Eqs. (Tb) -(10b) reveals that the operator 

_ determinant of the four simultaneous Eqs. (7a)-(10a) is of degree six. There- 

fore, the solutions (10d)-(7d) contain only six independent constants which can 

be used to satisfy the boundary conditions. The remaining eight can be ex- 
_ pressed in terms of the six independent constants. _ Such expressions, eight in 
"number, , are obtained by substituting the general solutions (10d)-(7d) into the 
original equations (7a)-(10a) and matching like terms. For convenience in 
carrying out these operations, the original transformed the 
A notation. For example, (8a) was rewritten as: 
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‘The general is now w complete. may be obtained 
ie inserting the boundary conditions for a specific problem into the oe 
solution. Some care is required for this step. Therefore, the constants will 
be evaluated for a particular problem to illustrate the method and give saahihe 4s, 
| 
fora popular structural system. (see Fig. 2). The ordinary simply supported © 7 
truss is a mixed boundary | value problem. ‘Three conditions are specified _ alg 
deformations and the other three are specified forces which “tell” ba he solution — 
that the been discontinued. They be written: q 
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_ These results are checked by calculating deflections for a | short i 


ia (63205052 


3) = ¢ 


al tool for the structural analyst has been illustrated in two cases. The reader 4 
will readily perceive that the two practical (but rather prosaic) examples — 
presented here do not fully « exploit the use of the method in solving | one- saeco 
dimensional lattices. Rather it is hoped that the two cases here, chosen from 
different structural classes, will serve as heuristic illustrations which might — 
interest others in working more interesting problems such as lattices \ which | 
_ are curved in space and have unusual basic wees or carry loads in a different 
ots, 
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ANALOG FOR EARTHQUAKE YIELD SPECTRA 


_An electrical analog has been designed to simulate a damped single degree 7 
of freedom, elasto-plastic structure subjected to an earthquake. Yield dis- 
| - placement | spectra of the 1940 El Centro earthquake are shown. It is em | 


ed that this approach be the basis | of a more relevant analysis. 9718 golts 


Earthquake resistant design of structures still largely on an cope 
cal basis. A ‘complete analysis of the situation presents great difficulties. : 
_ Buildings have been treated as simple, linear, damped, oscillatory systems of 4 if 
one or more degrees of freedom. Spectra have been produced by several in- i | 

vestigators, , showing the build-up of acceleration in such a system when sub- 

jected to actual earthquakes. These spectra show that in a large earthquake 
and with a reasonable damping factor, structures with periods in the approxi- 

y mate range of 0.2 to 0. 5 secs will experience accelerations of up tol g. eat. 

j Structures would certainly not stand these accelerations without yielding ; 
failing. _ The question arises as to whether the assumption of a linear analysis 
is tenable. It is observed that in large ‘quakes structures are often stressed 

r to the point of cracking of concrete or of mortar joints but remain otherwise 

suggested that in these cases the ‘structures worked through 


ld le. 
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A building coc code based on 1 the linear spectrum is is 
; because of the cost of designing for such large accelerations. It may well be | o | 
o that codes used in practice which design for accelerations much lower than 7 
those predicted on a linear basis already implicitly allow for yielding of the — 
: ee structure without explicitly stating so. It is proposed that we accept yield and 


base our code on a consideration of the safe extent and implications of this u 

with the preparation of yield displacement spectra 
ofa simple damped oscillator ibjected to the N-S component of the El Centro” 
earthquake, May 18, 1940—the worst earthquake so far recorded. ee 


_ The mechanical system considered is that of a mass and a stiffness - ba 
_member having an elasto-plastic restoring force. The restoring force is pro-— 
portional to displacement up to the yield point and then stays constant. Ifthe 
velocity now changes sign the force displacement characteristic decreases 
_ along a line parallel to the initial linear portion. The hysteresis loop is thus | 
a parallelogram. The system is also damped both in the elastic and inelastic = 
ranges. It is assumed that whatever | causes the damping continues to do so 
after yield takes place. When such a system is subjected to an impulsive ex- 
Citing force of the nature of an earthquake, the initial impulses set the oscil- — 
_ lator swinging in its elastic range. . At some stage an impulse from the earth-— 
quake may force it beyond the elastic range into the plastic range. It ed 


set of such a system is a function of period, damping, yield value, and the 
a DA 


An electrical “geslog has been set up to measure the maximum yield dis- | 


; "placement, t , the total energy dissipation and 1 final permanent set of this m me- 
chanical system. , The analog is basically | ‘the same as that used by Murphy, 2 
Bycroft and Harrison (1) for finding stresses in multi-storey structures. 

2 has been ‘modified in that biased diodes have been | placed in parallel across 
active condenser to ) produce the analogy of the mechanical hysteresis lo loop. 
The basic circuit is shown in Fig. 1. The condenser diode branch behaves as 

a condenser as long as the applied voltage is less than the biasing voltage. Ps 

‘If an attempt is made to raise the voltage above the biasing voltage current 
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x = x= (x~Y -Y) 


> 

where e is the algebraic sum of the charges which have flowed in the diode ~ 
_ branch of the circuit up to the time considered. : 
‘These two equations a are identical and the analogy follows. However, itis | 
> a necessary to work at much higher frequencies in the electrical circuit than in 
the mechanical circuit. In the system used the earthquake accelogram a(t) a 
was scanned p times as fast as it actually occurred and the resulting voltage } 
integrated. A current I(t), proportional to this — voltage, was fed into 


and (2) may now be rewritten to give -_ 

a(t) +4wrAx = 

= 


) > 


that 
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me This permits the analysis of a structure of varying frequency by changing - 
- only the scanning rate of the earthquake in the electrical circuit, = = 
i _ The displacement of the mass at any time is the sum of elastic and yield 


tracted from the maximum over-all displacement measured from 
a initial zero that occurs during the passage of the earthquake. It is the 
: maximum value of the algebraic sum of the yield displacements where 
- signs 0: of the yield displacements are taken to be the same as the veloci- 


Total Yield Excursion 

a ee is the sum of all th the yield displacements added | irrespective of 


tailed account of the driving and recording parts of the analog. ie 
? oa Briefly, the earthquake accelogram | is mounted on the outside of a 30 i 
diam. revolving drum. One side of the accelogram is blacked out and the a 4 
whole scanned by a light slit and photocell. The resulting voltage, proportion- ei 
al to o acceleration, is amplified and integrated. A. current I proportional t to g 4 
i ‘this voltage is fed from a pentode into the yield circuit and the voltage on the iS 
respective measuring condensers is photographed on a cathode ray oscillo- 
~ scope (C. R.O. ). The wheel is speeded up to 600 r.p.m. . and allowed to ) decay a. 
under an eddy current brake. The response of the circuit for each passage of 
the earthquake is photographed on the C.R.O. The film is moved forward each 


5 
erm of the wheel by a uni- -selector switch driven from : a timing mark on fee 


the drum. The slowing down of the drum changes the parameter p and permits a : 
an analysis of structures of different periods. _ The decay of the wheel is slow — 


avery detailed spectrum isobtained. 


= 


. A few modifications to the original analog have been made. A phase in- Brie 
_ version and switching stage has been incorporated so that for certain of the ie 


- the next presentation of the earthquake, the measuring side of the condenser a 
gg shorted to the high tension (H.- T.) by relays after each passage of the ro . 


cathode follower stage has been inserted between the measuring 


"and the so oscilloscope presents an impedance — 
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"ELECTRICAL ANALOG 


megohms. Ap photocell operated from a timing mark on the 
wheel triggers the ‘relays and the camera drive. As in the previous analog, 7 
minute marks appear on the film from the second beam of the C.R.O. The - 
distance between these minute marks on the film gives the period of wal | 
structure at the middle point ofthe marks, 2 
_ Fig. 2 shows the circuit for measuring the maximum yield auattesgent| 7 
Each alternate ee signal is reversed by the phase inverter and fed or . 
| 1 A D.C. current eats ps the plate of the pentode through the coil to H. *. 
and as the coil has a sensible resistance the D.C. voltage drop across itis _ 
balanced by a battery. Final balance is secured by the small variable balanc-— 


ing 1 resistance. This ensures that the diodes a are equally biased. es 
= coil of inductance, 24 henries 
= - damping resistor, 100 ie to 30 kohm 


= ae battery of 1 12 volts aa 
= resistor 0 - 500 | 

= of 0.0 003 to 0. uf 


measuring ‘condenser 1 16 uf 


= switch on measuring 


measuring condenser 


lar condenser being selected in any particular circumstances so that the aD, 
voltage drop across it is small compared to the biasing voltage. Its a, 
in the circuit i is such that it measures the total charge that has flowed in the 


diode branch and thus is a measure of the yield displacement at any y instant. i: 


Fig. 3 shows the arrangement for —- the total yield excursion and the 


_ The earthquake is fed in the s: same way — the time and the switching stage 
83 is used to switch alternately between the measuring condensers Cp. “di 
these condensers measure respectively the current through each diode they _ 
add up the total yield in one or other | direction. The switching stage adds them 
positively on the oscilloscope and gives the total yield distance irrespective of x 


sign. Subtraction of the yields in the two directions gives the permanent set. a i 7 


Switches S; and S2 discharge the pinata condensers after each passage of 


od, Before each run the resistance Rp ar and battery ep were adjusted so that the 
resistive drop in the inductive branch was balanced by ep. An infinitely stiff 
structure will yield when the base acceleration is equal to the ougivaient _—e 
ae acceleration Gg. . This may be used to set up | the gain n control. _ 
Fifty-nine volts were chosen to represent the maximum acceleration | of the 
El Centro earthquake i.e. 0. 32g. This infinitely stiff structure is simulated | Fis 
by rem removing the resonant condenser Cy. The diodes are then biased aad 5 
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is the required gain setting. The bias welne is then changed to 596/0. 32. 

volts to represent the particular yield value Grequired. = 
Calibration is effected by measuring the voltage across the inductance 


with the diodes and resonant condenser C,. Temoved. ‘This is photographed at 
the beginning ofeachrun, 
=a tt is necessary to determine the relation of voltage on a ‘measuring con- 
denser with displacement in the structure. Imagine the resonant 
C, to be very small, i.e., Co. This corresponds to a building being very 
rigid—the conditions of { calibration. 1 If Ko is the s | stiffness of this very rigid 
building and Ag the maximum ear thquake acceleration, aad maximum relative — 


is shown as the calibration height Z on the film. ons 


bigger measuring condenser Cm tl the displacement S would show as a te 


ioe All components : are accurate to within 1% and the over-all linearity of the ¢ 
amplifiers is within 3% from 40 cs/sec to 10,000 cs/sec. The behavior of the 
analog was checked by measuring the yield taking place when a single square 
wave is applied to the system and comparing this with the theoretical value. Be 
The results are shown in Fig. 4. The theoretical and experimental ratiosof 
yield displacement to elastic displacement : are shown as y /t and y/e e ye plotted — Sy 
against pu the ratio of pulse length to the period of the oscillator. Two values | oa 
of the ratio of height A to yield value G 


are shown in Fig. 5 for two structure periods are 


EM 4 
| — 
Aa 
at 
| 
it follows that a height y on the 
4T. 
vite 
ae 4 


7 lo WAVE CHECK OF a4 
. Figs. 6, 4, 8, and 9 show the maximum yield displacement taking place for — 
various yield values G damping coefficients d and periods T. The maximum 
elastic displacement X, corresponding to the particular yield value is also a 
shown. An addition of the two displacements gives the maximum over-all dis- 
placement from the initial 1 zero. Figs. 10 » il: and 12 give the total yield « ex- 
cursion and if this distance is multiplied by mG it gives the value of the Sa 
_ plastic energy dissipated in the structure. 2. The permanent sets are given in , 
"These results have several interesting features. ws It is, of course, to be ex- 
pected that increasing the damping and yield values will greatly decrease the 
_ yield displacements. If a building is designed to the common rule of with- a 
standing one tenth of its weight horizontally, safety factors will determine its 
_ yield value probably between G = 0.175 and G = 0.24. . Practicaldamping 
ratios A generally lie between 0.1 and 0. Z.: Figs. 8 and 9 show that the upper 
and lower limits of the maximum yield displacement occurring in this range = 
of parameters a are 1. 2 inches and 0. 2 these limits oc occur r at 


Be It may be seen from the figures that the period orresponding to the largest 
_ of the maximum yield displacements depends on the yield value. This period 
_ becomes shorter with increasing yield value. With the small yield value of 
_ 06 the largest yield occurs at about 2.5 seconds and is large. _ However, si 
_ ‘single degree of freedom structures almost invariably have periods shorter 
0.5 secs, and the yield displacement of up to ten inches is irrelevant. s¥ 


a The relevant maximum yield epee for Gs 0.06 will be less than two 


— | 

= 
= 

| 
generally have periods rather less than these, the significant maximum _ 
displacements will be less than one inch. A typical single degree of 
a _ freedom structure will have a about 0.3 secs and damping ratio 
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0. 06, 0. 12, 0. 75, and 0. 24 is respectively 1. 6, 0.6, 0.3, and 0.24 inches. In- 


tuitively, one feels that even a of inches in a wend of, 


fifteen feet is not excessive. 
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The te total yield excursion is seen to have the same general character as 


yields will have taken. place ‘aioe the. passage of the earthquake. 
in the case of G = 0.12 and a period greater than about two seconds, only ¢ one = 
_ yield occurs and this is shown by the fact that the maximum yield dis- 

_ placement, total yield excursion and permanent set are all equal. The | petal 
nent sets have the same and id in the region of 0. 3 secs period 
‘The extent to which a strectate may safely yield is a subject needing = 

_ tensive investigation. However, the other determining | factors such as per- 


a - missible movement of the structure and its permanent set uty more — 
Figs. 16 and 17 compare elastic and elasto-plastic designs. 
Two yielding designs of A= 0.10, G= 0.175 and A= 0.20, G= 0. 2, 
_ been chosen to represent the probable upper and lower limits of these parame- 
: ters. The figures compare the accelerations built up in the elastic case with | 


— 


the yield value acceleration Gg of the inelastic case and also the resulting dis- 
16 shows that a | logical ona linear basis would have 
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: 


to » stl for accelerations of up to 0.52 g. _ Fig. 17 shows that if yield is per- 
mitted the total movement is. about the would occur ina linear 


- Research is needed to determine permissible yielding of plastic hinges. : 
When this is decided an analysis of the type described may be used to design = At : : 
a a structure in a consistent manner » The paper at sie to show that many be 4 
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= TOTAL YIELD EXCURSION 
to the basis of 0.1 acceleration m must. yield in in large 


of ¢ earthquakes but that the consequences of this are not particularly serious. a4 


as The total movement is often less than if the structures had been designed © 


= up to 0.52 g. The extent of the yield displacements is generally less _ 
_ than an inch, and an inch yield movement of one floor to another represents ~ 


4q an extremely & small plastic strain in the plastic en. The permanent sets 


— «3. | 
— 

— 


ay) 


in Inches 


no 


Permanent Set 


= 


> 


4 | 
| 
= 
a 
4 
| 
| 


1950 


Inche 


ermanent 


P 


7. 


aiding Period Secon 


AMANENT SET. 


=O-175, 


plastic hinges are available. ‘This applies particularly to steel framed — an 
: ee However, if the yield value was lowered too much the res resulting 
increased displacements would probably be the limiting factor. 7 
De ai Other building types such as bearing wall buildings | or buildings involving © 
pre-stressing do not fall into this | category. it would appear logical that these | 
types should be for accelerations of up to half that of gravity. 
cy. Dr. . B. W. Faulkner(2) ‘kindly eo out to the authors | one result of allow- 


7. : = likely that design coefficients could be safely lowered in buildings where 
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superimposed on the permanent oat caused by the asymmetry of the earth- il 

— This feature may increase or decrease the fundamental permanent set. _ 
‘The yield spectra of more earthquakes must be determined before definite eee 

conclusions are reached. It is likely though that they will not differ substanti- _ 

ally from the one earthquake | considered here. 7 

_ The authors hope to analyse more earthquakes on the same basis and to 


extend the the work to yielding of malti- -storey building frames. 


‘: “Electrical Analog for Shear Stresses in a Multi-storied Building”, by 
M. J. Murphy, G. N. Bycroft, L. W. Harrison, Proceedings of the World « 

Conference on Earthquake e Engineering, Berkeley, June 1 
2. De, Pauthner, Design Bagincer, of Weeks, ‘Wellington, 
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THE NEUTRAL AXIS IN PLASTIC BENDING OF BEAMS@ ~ 


Aris 


7 


| 


7 The motion of the neutral. axis in ee pure | bending of beams in in 
plasticity is discussed. Both ergs plastic beams and beams with linear 
strain hardening are considered. A method is given by means of which it is a 
possible to o correlate in a simple manner the motion of the neutral axis with Y 
} the motion of the elastic-plastic boundaries. Formulas are found for the 
motion of the neutral axis with increasing curvature. Conditions are es-— 
tablished for reversal of direction of motion of neutral axis, 
perropucTion 

purpose of this paper is to discuss the motion neutral axis with 
‘increasing curvature in plasticity, and the relation between the motion of the - 
neutral axis and the corresponding motion of the elastic-plastic boundaries. | 
This discussion will be restricted to symmetrical pure bending of beams. e... 


Specifically, it is assumed that the cross section has only one axis of symme-__ 
try which lies in the plane of the external forces, and it is also assumed ‘that a fs 
the stress-strain curves for tension and for compression are identical and 
that they have the form shown in Fig. 1. On the basis of the elementary theory _ 
of bending of beams in plasticity(1) a simple graphical method will be outlined 
by means of which it is possible to find the elastic-plastic ‘boundaries corre-_ 
sponding to a given position of the neutral axis. In the last part of the paper | 
certain conclusions will be drawn concerning the motion of the neutral axis — 
and the motion of the elastic- -plastic boundaries for increasing value of the 
curvature, and formulas will be found expressing the rate of motion of the 
Note: Discussion open until March 1, 1960. To extend the closing date one month, a ; 
written request must be filed with the Executive Secretary, ASCE. Paper 2198 is — 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedin » 


of the American Society of Civil Engineers, Vol. 85, No. _ EM 4, Ontober, 1969. 
. This investigation is pert of a researc program sponsored by the National 
1. 


. Associate Prof. of Civ. Eng., Yale University, New Haven, , Conn. ‘. 
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axis. Also conditions will be given for the Teversal of the direction of motion 
In Figs. 2a and 2b a general cross section has a plastic region on on only | one 
- side of the neutral axis. In Wig. 2c the cross section has plastic regions on 


7 neutral axis as a function of ‘the curvature and of the position of the ‘neutral 


stress distributions are the same Figs. 2a to 2c. Because of 

; symmetry the neutral axis is always normal to the axis of symmetry. pee” ail 

_ For a small bending moment the entire cross section is elastic and the A. 
neutral axis passes through the centroid of the cross section. The fibers hy 

which lie at the greatest distance from the neutral axis will always have the 


ning from these fibers. These fibers will be either the lowest or the srming begin~ 
fibers of the cross section taking the axis of symmetry of the cross section to 
be vertical. if the highest and lowest fibers have the same distances from the 
af - elastic neutral axis then they will all become plastic simultaneously. If the ; 
a ; highest and lowest fibers are at different distances from the elastic neutral 


— axis s plasticity v will start at the fibers which are situated at the , greatest _ sit 


2 is only a schematic representation of different cases of discussion; it 
does not represent r r >sults of actual calculations. 
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| distance from the elastic neutral axis. - - what follows it will be assumed that : 
the fibers situated at the greatest distance are the highest fibers; hence the _ 
cases given in Figs. 2a and 2c will be considered. Obviously, a rotation of a . 
the cross section by 180° produces the case 2b from the case 2a. Thus by re- 
stricting our considerations to Figs. 2a and 2c we do not restrict the validity 
of the argument to these two cases only 


slope of the line OA in Fig. 1. This the modulus: of elasticity, 
* slope of the line AB in Fig. 1. This is the hagas af of linear 


we 
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= distances of the neutral axis from the lowest and highest fibers a 


respectively. the position of the neutral axis changes with» 
increasing curvature (once plasticity has started) it follows “4 


= total height of cross section 
yiel strain 


oa strain of the highest fib eel (considered positive) © 
= curvature 
= distances of any fiber from the lowest: and highest fibers re- 


dl spectively. Referring toa given fiber these di distances are 


ny, ny! = distances of lower and upper -plastic boundary from 
lowest and highest fiber respectiveiy. These distances are vari- 
able; they are respectiv ely zero wnen the lower and upper = 
Bin), Bir!) = width of the cross section at distance n or rn fro low = 
cur 


A _—__ = total area of cross section. 
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‘Thee expressions for r Ag(n) AS ) have been calculated some cross 
_ sections and are given in reference (1). For example, fora rectangular cross 
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plasticy Fig. 3 the graphs for the six ‘functions B(n'), A(n), A'( Aa(n), 
) are shown for a general cross section. slope of curve at 
lowest | 7 = @, with respect to the vertical is equal to the ordinate A(a) and to the area 
od of the cross section between 7 = 0 and 7 = a. Similarly, the slope of the _ 

7 - curve A9(n') at n' = a', with respect to the vertical is equal to ~ ordinate 
a A'(a') and to the area of the cross section between n' = 0 and n' = at. Re- 


= mark also that as (or n') increases, the value of A(7) (or of Ad( 
creases also and thus also the slope of (or of Ag(n')) increases. For 
a= 0 the area is zero; hence the curve Ag(n) is tangent to the vertical at 
= =0. Similarly, the curve Ag(n' ) is tangent to the vertical at = 0. 

aaa a of the neutral axis will be e given | by th the expression 7 


a As it is seen from Fig. 2c, when plastic regions exist on both sides of. the = 


neutral axis, the distance of the upper elastic- from the 


cross 
r cross 


*Expression (1) follows from formulas (4-19) in reference (1) by ie pit Pov 
Np = Na and noting that in the present case there exist only 


bers 
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In Fig. 4 the lines s Ao(n) and AS(r! ) are e drawn for a specific cross section. 

_ These two lines meet at Ny' where Ag(7) = Ag(n'). Note that when no ar 
exist—that is is, in the (1) becomes 


(h )= = ( h’) 


‘then = Ag(nj) = 0. Hence the elastic position NNj of the q 
neutral axis passes through the Nj of of the two lines Ag(n) 
is drawn in Fig. 4. This line starts a the inwouetinn N between nelastic 
neutral axis and vertical base 00' and is drawn toward the boundary which is , 
at the shortest distance from NNj. . He Hence, in Fig. » 4 this line NM is drawn , 
toward the lower boundary. —iItis seen en that 
ab) = —— a 6") 


‘because (a'bt) = - ‘Therefore ab 2 a'bt, with the equality sign 
being ve valid when Ej - = 0, that is for a perfectly | plastic ‘material. - On the — 

hand E; = E fora purely elastic material; hence for an elastic material a 
a ab = ~ in which case se the curve MN coincides with the line NN]. ee. =: _ 
‘The first a appearance ofa plastic region will be at : the boundary of the cross— 
_ section which is at the greatest distance from the elastic position of the 
neutral axis; in Fig. 4 this is the upper | boundary. _ When a plastic region ex- 
_ists om, one side, for example above the neutral axis, then in Eq. (1) we 


ad Let nn! be a neutral axis in Fig. oe ‘ai the left hand member in Eq. (5) 
z is ‘equal to length ns and the right hand meee will be found by drawing a 
vertical line from 8 until the / A» 
exists only on one side of the cross section the corresponding values of h' and 
Ta "1 will be determined by drawing vertical lines and finding the intersections 
oar of these vertical lines with the lines NM and nition The intersection of the a 
vertic: 
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_ We remark that length O'p = 7, length O'n = h', and length On = h 
increases it is seen that h incre reases also but at a ‘slower rate because the = mt 
slope | of NM with respect to the vertical is greater than the corresponding 2 
slope | of O'N]. It appears as if the neutral axis tries to avoid the plastic rie 
regions as this plastic region spreads downward but that it does not succeed » a 
to keep the distance from the elastic plastic ‘boundary « constant. From formu-_ 
la (4) we also see that 
‘which aids that the distance between the elastic-plastic t boundary and the aa 
neutral axis is inversely proportional to the curvature k. Expression (6) can te 
be used for finding the curvature k which corresponds to a given position of. 
pa The plastic region will be restricted to the upper side as long as length 
sr ~ length st. ‘Suppose that for the position mq of the ‘neutral axis it sl jot 
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= will be two plastic regions, one above the neutral axis and another below 
aye Consider now a position cd of the neutral axis below position mq, 
in Fig. 4. Suppose also that line eg is the elastic-plastic boundary above the a | 
neutral a axis and line hi is the elastic- -plastic boundary below the neutral axis. 


Thus we have 
Consider now the line kg parallel to 00' in Fig. 4 and the extensions of line 
hi to k and of line cd to u. . Obviously ku = ug and ik = : ed. ‘Suppose that kg and 
ik are the two normal sides of a right angle triangle with the right angle (gkh). 
_ Let this triangle be moved down and to the left in a translatory motion so that 
the length ik will remain constant and let this motion be continued until the a 
length ik takes the position 0k'. Then the vertical side of the triangle will 
_ pass by d, so that the new position of the triangle is Ok'd. ——- 
i Suppose now that it is asked to find the positions eg and hi of the itis 
plastic boundaries corresponding to the position cd of the neutral axis. ‘The r 
5 translatory motion mentioned above is reversed so that now a right angle tri- 
age is used, , the vertical side of which passes by | d while its horizontal side 


now moved in a translatory motion so that the length Ok' remains constant — ' 
until at some position (ikg) there is ku = ug. At this position hk and eg are fae 
the two elastic plastic boundaries and Eq. (6) gives the corresponding curva- 
ture. _ Fig. 5 shows the procedure and in Fig. 6a complete determination of a 
the: elastic plastic ¢ boundaries is given for several po positions of the neutral axis 


€ 


= a section and a given material. 
Discussion 


"any point s ie always greater than the slope of the line O'Nj at the poetic 
a point r. Hence when only one plastic region exists the neutral axis moves 
always downward. When the neutral axis intersects the MN line at q the lower 
elastic plastic boundary is at 0. But the slope of ON] is zero at 0; hence it is 
_ smaller than the slope of O'N] at l. Therefore, when the second plastic region | 
_ first appears the neutral axis continues to move upward. 
_ The motion of the neutral axis as the curvature increases can be related 
to the motion of the elastic-plastic boundaries as follows. From Eq. (4) we | 
*Here q is the intersection of line if with line NM, and line mq is 2 a neutral eS 
axis—not shown in the passes by the point q. 
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(10) gives the > rate of motion in terms of the rates 

of motion of the elastic-plastic boundaries. It is seen that for the neutral axis: 

to remain stationary the the ‘rates of motion of the two boundaries must be equal. "7 
It is seen also that the neutral axis ‘moves toward the elastic-plastic boundary ¥ o 

a Eq. (10) is not sufficient to give the rate of motion nof the neutral axis be- ie 
cause the quantities dn/dk 2 and dnj/dk are not known. Curves A(7) an and A(nt) 
can now be ow be used rate motion ¢ of tl nail neutral axis. Eq. (a) 
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“while the right hand of Eq. (11), with the help of (9), be 


— 


Recalling that dA» 9/dh is area cross section below the neutral 
axis and dA3/dh' is the area of the cross section above the neutral axis it is _ 
concluded that the sum of these two quantities which appears in the denomi- - 


nator of is to the area of the cross section. 


#8 is recalled that dA2/d n1 is is the lowe lower pl: plastic area of the cross section 
and that ‘dAo/dn: nj is s the upper plastic area. As EA/(E - - Ej) is always larger > 
or equal to A it follows that the denominator of the right hand member of 
Eq. (13) is always negative. Therefore dh/dk has always a sign opposite to " 
that of the numerator of the right hand member of Eq. (13). It follows a. 
fore that the neutral axis moves upward, remains stationary, or moves dowa- - 
ward when dA3/dnj{ is smaller, equal, or larger than dAo/d "1, , respectively, — 
that is, when the upper plastic area is smaller, equal, or larger than the: lower uh 
Plastic area respectively, 
When plasticity first appears the upper plastic area is the only 


one; hence, the neutral axis moves downward. This downward motion of the 
neutral axis continues even after the lower plastic area appears. As the eld j 
plastic ; region increases in size it is possible that this region may become | 
equal in size to the upper plastic region in which case the neutral axis stops 
‘moving. When the lower plastic region becomes larger than the upper rate 


tf the ie slope i is greater than the slope at. g the axis 
upward. 


— 

i — 
q 
| 
Consider also the positions (cd), (eg), and (hi) of the neutral axis and of the 
os _ two corresponding elastic-plastic boundaries in Fig. 4. For the neutral axis _ 4 Ge aa 
ia remain stationary as the curvature increases it is necessary that the slopes 
_(with respect to the vertica]) at i of Ao(n) and at 2 of A4(n') be equal. Ifthe | 


Eq. (a3) together . with | Eqs. (8) are are sufficient for determining the rate of 
4 ‘motion of the neutral axis. Indeed Ag(n). and amy are known functions and 


, & Thus, at each value of x and h we find from (14) the value of an/ax, pee which 
me: the new set of values h + dh and x + dx can be determined. Thus, a step by a 
step integration « of the differential Eq. (14) can be performed. ‘The graphical 
: "construction given in this paper seems, however, much simpler to use than > 
hs Strain hardening influences the motion because of the factor E/(E-E,) in 
_ Eq. (1). For Ey = 0 there is no strain hardening. For E = Ej there is linear 
"elasticity. As strain | hardening increases E-Ej, b becomes smaller and the al 
middle line NM in Fig. 4 tends toward the horizontal position. Hence as strall 
; hardening increases the motion of the neutral axis tends to slow down for the 
; same amount of motion of th the elastic plastic boundaries. It is also seen that 


two materials with the same ratio E/ (E- — have the ain a be- 


value of ane E,). Compared to the well known method of Saint-Venant de- 
ale in ea (2) the method presented in this paper has the advantage | 
eo. in order to find the position of the elastic plastic boundaries correspond- 
ing to a given neutral axis now only a comparison of lengths is ene? § | 
_ in Saint-Venant’s method a comparison of areas is necessary. i 2 a 
Ino order to find the relationship between the motion of the neutral axis and 
the motion of the elastic ~plastic boundaries it is necessary to repeat a given | 
procedure a number of times. Such a repetitive procedure is much easier 
 seemgliiet by means 3 of the present method than by means s of a method in 


Phillips, “Introduction to Plasticity”, The Ronald Pp 


, vol. Il, Third Edition, D. Van 
4 Nostrand Company, York, 1956, 366. aed 
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LATERALLY LOADED THIN FLAT PLATES 
illiam A. Bradley, ! M. ASCE 


SUMMARY 
In the work reported in this paper, certain. laterally loaded thin plates were 
‘studied experimentally, using the Moire method, and analytically, by finite ay 
‘differences - The procedures followed in the study are summarized, the re- . i 
_ sults are presented, and possible sources of error in the Moire technique are 


discussed. In general, the agreement of deflections and moments as found by | 
q the two methods is good, maximum variations being about ten percent. bs 
 ‘Itis concluded that the Moire method, with possible improvements in tech- 


nique, offers a useful tool, along with numerical methods, for the future study 


es the formal analysis for - moments and deflections in laterally loaded tl 
plates, the formulation of the boundary concitions and the solution ¢ of the p =. 
tial differential equations, except for certain of the : ‘simpler cases, become — 


difficult. of the alternate procedures which may then be to 


— 


< ‘be used to determine expressions for moments and deflections. — 
He Pe Numerical solutions of the differential equations —in which the differen- 7 
tial equation can be approximated at a ame number of points in the 


i 
plate by finite difference expressions. ‘This procedure leads to the 
solution of systems of linear algebraic equations—a routine which can 
ke cae by digital « computer or byt the e use | of | the desk calculator. — 
Note: Discussion open until March 1 1, 1960. To extend the mn date one mau. a 
- written request must be filed with the Executive Secretary, ASCE. Paper 2199 is 
wc part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings — 
th the American Society of Civil Engineers, Vol. 85, » No. EM 4, October, 1959. 
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‘methods and the results compared. The experimental procedure used was 
the Moire method, a photographic interference method first described by ve 4 
Ligtenberg, (1) and later applied by Bradley(2) in the study of clamped pl plates 
. To study the application c of the Motré tasthod to plates : 
se those with simple supports, in order to to determine its enmeel I 


determine the deflections the of moments in certain 
_ plates by both the Moire method and by finite differences. The com 
parison of the results by the two methods serves as a check on both. 4 


_ Some of the results obtained in this study have been reported by B. Raju - 


with the shape of the ‘deflected plate being given the oft the plate 


in which D, the plate constant, is to bei 
_——s- The derivation of these formulae assumes that no normal stress acts in 
; the middle plane of the plate; i.e., no membrane stresses exist. ‘For small 
a _ deflections, this assumption is acceptable, but if the deflection is of the order 
‘of half the plate (Gicknese, these membrane stresses cannot be be neglected. 
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‘a 


find w = = f(x,y) such that equation (6) and the boundary conditions of the plate a 
are satisfied. Knowing the function, w, throughout the plate, equations (1) _ a 
through (5) can then be used to calculate the moments and shears. es uh 
If the solution of equation (6) was unobtainable, it would also be possible to — 
determine My and My at at any point in the plate if the curvatures, 


nc 
aya _ the plate constant, and Poisson’s ratio were | known. These curvatures: 
were experimentally determined in this by ti the Moiré method, which 
gives the slopes of the plate at every point in any desired direction; from cS : 
these slopes, it is then possible to determine the curvatures. J Another experi- 
mental method, in which curvatures are directly measured, is applied “ere 7 
mercially by the Presan of slabs. 


_ The Moire method is explained in detail in both Aone lis 1 and 2, and is 
summarized in reference 3. The essential items required are those shown a 
in Fig. 2—the model having a reflecting surface, a finely lined grid of alter-— 

nate dark and light lines of equal width, and the camera, located behind a a 


small hole in the grid; the test setup is shown also in Figs. 3a and 3b. An = | 


exposure is first made on the negative with the plate model unloaded, or v with bs ll 
a small initial load; “next, the final load is applied, and a second exposure is por 
made on the same negative. The interference of the reflections of the grid = ‘ 
produces patterns as as shown in n Fig. 4. These patterns give the : slopes at all 
points in the plate model; for example, in Fig. 4a, which was taken with the 
grid lines in a vertical position, the fringes are those for slopes of the loaded 
plate in the horizontal, or or x, direction. Since the top edge is clamped, th the 


slope along this edge i is z zero, and the fringe is of zero order; a at each point — i 


along the next in x direction be punity 
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‘Fig. side view of test set - 


"multiplied bya a constant which is shown in reference 2 to at e equal to 5. a — 
: /4 where d is the spacing of the grid lines and s is the Gatanee from the model 
to the grid. Thus, in Fig. 4a, the maximum slope of the model in the el 
 Ganstion is (6- 1/2 ( es ) at the lower right hand corner, one in Fig. 4b, the 


maximum slope of the model in the y direction is (13- -1/2)(--) at the lower 
After the slopes in a given direction are : found, it is then p possible to plot 
curves of the slope along a given line, as shown in Fig. 5. Then, the anpe of 
- VS X curve at a given point, Pp, in the plate will give the ——— value 


, the twist; from th this value twisting moment, 
is also to plot values of of Ja along a given and 
> ea the : slope of this curve will give the derivative sndieiaeh in 1 equation (4) or - (5) 
for computation of the shear. It should be noted, however, that this has cer- 


tain practical limitations, since graphical produces errors and 
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a The deflections may be aie by integration of the area under the slope 
“ _ The relationships between the slopes, curvatures, and moments in various 
al directions at a point in a deflected plate may, in some cases, be useful. Using — 
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concentrated load at lower right hand corner. 
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series of fringe photos for a uniformly - -loaded plate is in Fig. 


te fringe patterns for become somewhat ‘obscure along the ‘side edges. 


Rweven, te nee fringes are well defined along these edges; thus, the 


“curvatures in the diagonal could be used along with (whi ch, 


this this case, i is the x curvature ( ) by ation of 


“equation (11) or (13). The use of the diagonal | relations is also helpful as a 


‘The solutions of differential equations - —both ordinary and partial—can often 
“be: approximated by the use of finite differences. This procedure is discussed = 
in detail by Salvadori(5) and others, and its application to plate problems is 
reviewed in reference 2; in ‘reference 2, also, the difference operators for = "ha 
points near various types of boundaries are summarized. The systemsof 
algebraic | linear equations which result from the application of the difference | &§ 
operators were solved largely by the Michigan State University digital com- ; 
puter—the MISTIC; a few of the simpler systems were solved by the Gauss- > — 
Doolittle reduction scheme, (6) using a desk calculator. Si i 
an For most of the cases given here, difference solutions have been anectet et 
out with two different grid spacings, and Richardson’s were used to = 


produce 


My = — + cos 26+ Myy sin2@ 
a 
as 
— 
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. 7. Moiré fringes for clamped, uniformly loaded square 


plate with cutout. Plate mode, 9"' x and 


? 
Experimental ¢ Considerations 


-arrangements—is Pretond in Fig. 3. “The setup is described in detail in 
references 2 and 3, but certain f features are 


* _ Both plastic and metal sheet material h has been used for the models. at 
vantages of the plastic 1 ‘materials’ 


a Smaller stiffness--thus making it possible to use thicker plates with — 
small loads and to allow greater deflections, with the corresponding 
finer fringe patterns, without deflecting the plate to the extent that 
8 large membrane stresses occur. The relatively small modulus of 
Poa elasticity of most plastics also makes it possible to obtain clamped — 
edge conditions more easily. Also, since small loads can be applied, 
_ the stiffness of the supporting not ‘ates to be be as great 


to of the frame. 
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be Surface smoothness—many plastic are poured on plate glass, 
i _ thus giving a very smooth surface, which is necessary for faithful = ae 
7 reflection of the grid. Ih rolled metal sheets, very slight roughness _ 


besa in the roller surface is transferred to the sheet. ae al 


-Machinability — plastics may be and d accurately ‘and 

Thickness uniformity —thie is ‘helpful, since any variation in ‘thick- 

ness will change the distribution of moments in the plate. Also, as 
‘shown in Equations (1) and (5), the > moments and shears at a point in the ~ 
plate are dependent on D, , which in turn varies as the cube of the thick- 7 


ness. This makes it necessary to accurately measure the thickness of 


sass as constant for metals, ‘but E for plastics varies materially with both 7 
temperature and humidity and somewhat from sheet to sheet. This 
: makes = necessary to calibrate the model each time a test is run. In 


elastic material. In using plastics, the creep tendencies of the a = 


torte must be considered in planning the experimental routine. 


_In view of the advantages of the metals, which it was felt would result in ea 
more accurate results, considerable effort was directed toward the use of eo 


to the eye ‘still peasensed enough grain to cause distortion of the grid reflec- 
tion. This, along with the consideration of the other listed factors, led to the 
adoption of the plastic for most tests. The material used was Perspex, .. 
methyl methacrylate; the transparent material can be used, but requires mir- 
oring of the reflecting surface; to avoid this added step, the black opaque 
sheet was used, and it yielded ‘good fringe patterns. All results given here 
were obtained using the Perspex model unless specifically stated otherwise. | 


| 
The number of fringes appearing ing for a “a given — plate depends on the 
‘Spacing of the alternate black and white lines of the . grid. In most of the a oe 
so far done, the line width is 0.05 inches, giving a d of .10"'. Since the dis- 

tance from screen to model, s, for the setup used is 25"", the difference in 

slope between two adjacent fringes is d/2s = 0.002. An alternate grid was 

used, with d = .05"', but it was found that the definition between a wes a 

Support of the Model 34 


In this study, the plates have been or 


With free edges. In all cases, the clear we 9. 05 


al 
was placed between the two 3/4" thick steel plates and the 8 the model - te 
together with 16" bolts. This assembly was mounted into the frame 
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_ For the simply supported edge, the fixture was as shown in Fig. % the 

, ‘model was held vertically between the knife edges by two fine threads. — The 


bolts were tightened only enough to lightly draw the assembly together. On > 


loads were applied to the means of a lever system, 
which appears at the rear ofthe modelin Fig.3. 
. > ; The un uniform loads were applied through air cells placed against the back 
7 surface of the model; one of these cells ready to be put into place is shown | 
in Fig. 8. The load pressure was read, in inches | of ae from the the a 


on frame as shown in Fig. 3. 


‘Fig 8. Clamped edge fixture; below is the air cell used for 
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Calibration of the Material ws sath, ory 


Since the ‘modulus of elasticity ity of Perspex varies with both ‘the tempera- 
ture and the humidity, ‘it is necessary tod determine it during each test. This © 
can be done by measuring the deflection, slope, or or curvature of a plate. of ‘_ 
known dimension under a known load, and, using the theoretical relation be- 
tween these quantities, calculating D. - In this ‘study, the calibration plate tion 7 
cut from the same sheet of material as the model, thus avoiding | the variation - 
in| material which may occur between different batches. <r 

_Two simple cases of plates which may be used for are dis- 


= a) Square plate supported at three corners and with temrensesnin load at an i 

the fourth. (Fig. 10 Pay 

this plate(7) the curvature in the x direction is constant throughout, 


= 


3 


a 


but, since between two su fringes, changes 


a P (2s) (spacing of 


Thus, knowing d and s (which are constants), the applied 1 load, and Poisson’s 's 7 8 


ratio, itis possible to measure the spacing from the Moire as 
in Fig. 10b, and i compute D. 
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This procedure that the calibration plate be mounted in 


the test frame and that photos be taken ie each test. Toa sie this, th the 


rat 
-* ‘For a uniformly loaded ‘iautid circular r plate of radius, a, ‘the center de- 


Thus, after loading the plate and measuring deflection, the 
con stant can be computed. . This measurement was made using the setup 
Asa partial check on this ouieuiiaie, ‘the Moiré fringes were obtained for 
‘this plate, as shown in Fig. 12. . From these fringes, the center deflection was 
computed as 0.0111 inches, and using equation (21), De = 82.3 Ib. -in. The 
_ curvature was also found, and d using equation (2), the center moment, My, , was 
“calculated as is 0. 549 in. -Ibs. p per inch, using Poisson’ s ratio as 1/3. This com- 
>» complete the eae the dial ga ones ia was then read for the. same 
a load, with readings of 0.0113" and 0.0110", or an average of 0.01115"". 
ol For all cases, the calibration plate v was chosen from the large Perspex: 


sheets so that the thickness variation w would be as _— as — 


sine Since the model materials used are plastics and subject to 
the accuracy of the results probably depends somewhat on the time intervals — 


during the test. . The procedure followed in this study was essentially that . 


photos to allow recovery ‘of the material. 
allt fringe pictures were taken using Royal | en wisi film 


Me 
d Reduction of the Data 
7 = After colangementa © were made of the fringe patterns for slopes, curves of 
the slopes were plotted in each desired direction by locating the half-fringes, 
since these are more sharply defined than the integral fringes. . From these, 
_ the curvatures were determined as indicated inthe Section “The Moire 
- Method” and Fig. 5. The slope curves were drawn on 20 square-to-the- -inch 
cross section paper, and the slope of the slope curve at a point was taken as 
the slope of a line intersecting the curve two small divisions each side of the 
Bsc speed the curvature determination, a special protractor-like device 
_ was constructed which was calibrated so that curvature could be read directly. 
The value of the plate constant, D, was found for a particular point from > 


me 
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Al Ae 11. Set up for calibration of material. The box appear- 
in background houses an air reservoir used to 


cushion im act effect of load lication. * 


pid 


Sen 
‘Fig. 12. Moiré pattern for loaded clamped dcalibra- 
<2 
te was taken as the av erage thickness of the | calibration plate, and 


"was the thickness of the model at the point. 
The moments were computed from equations s (1), (2) and (3), using 


Fetesad? s ratio as one-third for Perspex, and equation (7) was used for « oo a 


termining deflections, | with the integrations being carried out by totaling oe 
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mare, Uniformly Loaded, Clamped Plate 
The fringe patterns for this plate are shown in Fig. 13. Figs. 18a and 13b 


show the slope fringes for a load of 0.2167 psi.; from enlargements of i 
3 these two photos, the moments were computed at intervals of 0.1la, a being the 
: 7 side dimension of the plate; deflections along the center lines were also coms 

: = In Fig. . 14, the moment values are given, as computed by the Moire 
method and also by finite differences, using a grid spacing of \ = a/10. In the 
case of the Moire moments, the four values given are those obtained at ‘sym- | q 
- metrical points in the plate and are presented to § give an idea of the range of E 
_ values. The difference results for = a/10 are slightly lower than the actual, § 
or extrapolated, values, as shown in Fig. 15, where the moments along a” 
-center-line are plotted; this figure shows values for A= a/8 and A = a/10, the 
3 extrapolated value, and the Moiré value. Fig. 15 also shows a comparison of 
deflections along a center- line. Ih Figs 13¢ and ‘13d, the fringe patterns for 


ae The curves of Figs. 14 and 15 show that the Moire and difference results 
, 7 are in very close agreement; it might also be noted that the maximum deflec- 
i tion is 0.0233"" for the smaller load, or about 0.182 times the plate reas 
= When the load was doubled, the deflection was 0. 0450"’, not quite double the 
_ previous deflection; for this load, the maximum deflection was 0.350 times 
i the plate thickness. fe For the larger load, a part of the load is apparently re- 
sisted by membrane action of the plate. ay: 
“Square, Uniformly Looted, Simply Supported Plate 


Fig. 16 shows the fringe patterns for this case, and Fig. 17 compares the 


 Moiré method and by finite ‘differences. The complete. distribution of My 
_ My and Mxy is given in reference 3, aie with the finite difference results - 
grid spacings of both a/8 and a/16. 


_ One of the purposes of this study was to determine whether simple support 
conditions could be obtained using the knife edge arrangement shown in Fig. 9. 
2 The boundary condition at the : simply supported edge leads to the requirement 
that the curvature be zero in the direction of the edge and also perpendicular 
to it. For the first requirement to be met, it is necessary that the slope in 
the direction of the edge be zero; looking at Fig. 16, it is seen that this is 
"particularly the case. To meet the requirement that curvature be zero per- 
_ pendicular to an edge, it is necessary that the slope perpendicular to the edge 
be constant; that this is true is indicated by the fact that the fringes (lines of 
constant slope) are indeed ndeed perpendicular to the edges at the intersection a 
The moment results obtained by the Moire method were, in general, five i) 
ten percent less than those obtained by differences (which, for this case, | ial 
agree very closely with theoretical values). The Moire deflection results _ 
were also about ten percent below the exact values. . Most of this beglanats. 
_ can be accounted for by considering the membrane | action of the plate. — 
a plate is deflected | Dy loading, a soutien the load is 
3 sisted by the flexure of the plate and a ine of the membrane stresses in 


— 

q 

— 
Pa 

— ii 
— 

— 4 


i a 


450 


plate, ‘Plate ‘Square; average thickness = 1281". 
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"uniformly loaded clamped plate. 
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the of small deflections, the effect of membrane stresset 

is usually neglected and the deflection is taken to vary linearly with the load, 

as shown in Fig. 18. If, however, the membrane stresses are considered, as 

» discussed by | Timoshenko(7) and by Prescott, (8) the relation between load and 

This relation for the plate of Fig. 16 was com 

- puted, and is plotted in Fig. ‘18, along with the approximate percentage of load 
‘taken bythe membrane action. = = = 

Ih the procedures followed, an initial load was always put on the plate be- 

the initial—or “unloaded condition” —exposure was made. In this case, 

the initial load was ‘equivalent to 2 inches of water, or 0.0722 psi., and the : 


- final load was equivalent to 4.5 inches of water. For a load of 4. 5" of water, § 


_ the calculated deflection, considering membrane stresses, is 0.041"" and for 
load, 0. ‘difference is 0. 0216", or 34.6 X 


‘model deflection for | this increment of load was 0.0230' 


- Prescott’s formula, used in calculating the deflections of he assumes ~ 
_ the model to be held, allowing rotations but no edge displacements in the 
plane of the plate, whereas the support used probably did allow some such 
; movement. It would be | expected, | then, that the model deflections would be 
_ between the value predicted by small deflection plate theory and the modified 
theory in which membrane stresses are considered; this was the case. 
a i a Inthe clamped plate of the Section, “Square, Uniformly Loaded, Clamped 
ys 3 Plate, ” smaller initial load (1"" of water) was used, since it was not necessat 
_ to seat the model against the supports; the increment of loading was 6" of | 
water (0. 2167 psi.) and the deflection due to this increment was 0.0233" or 
about « one-sixth of the plate thickness; the maximum deflection due to the 


total load was about 0. 027", or 0. 21 times the thickness. _ The greater deflec- 
tion, 
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7 a Simply supported plates with both rectangular and square cutouts and sub 
is ject to distributed loadings have been studied; the fringe patterns for one of 
these cutouts are shown in Fig. 19, and the comparisons of Moiré and dif- — 

- ference results are given in Fig. 20. These results were calculated using — 

~ the fringes obtained with q = 1.5" of water. Moment curves and deflection 


; values for other Cutouts, obtained both by the difference and the Moiré = I 


indicated that the effect of a cutout in a simply supported plate extended a 
distance from the cutout about to the 
End Supported Rectangular Plates oh 


= results obtained by the two methods were about 10 percent. "These studie 


with concentrated center load are summarized below; in both cases, the side 


edges were fr 


method, using | oqeare grids with a = and. and a = 12. 


= the values of deflections and moments : as determined by the ditferenc 


Moiré fringes. for simply supported plate with un uni- 
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Fringe patterns for plate, ‘clamped at 


| deflections along the centerline of the plate is given in ‘Fig. 25; and, the finite’ 
- difference moments and deflections for A = a/8 and for A = a/12 are given in | 
‘Fig. 26. Note that at the section midway between supports, the deflections at 
2 edge of the plate are greater than the deflections under the load. =~ ts 


this a square model of steel was used. The thickness 
was 0.038"" and Poisson’s ratio and the modulus of elasticity were determined 
by test to be 0.243 and 29 x 106 psi., respectively. The surface was hand- — 
“polished using al aluminum oxide powder; that this surface was no not of as — 
- quality as that on the plastic models is indicated by the fringe patterns in 
The initial the 2 pounds and the final 8 pounds, 
with the incremental deflection between these loads being 0.0398", or slightly 
greater than the plate thickness. . As expected for this large deflection, the 
actual value is less than calculated by finite differences using the small de . 
- flection theory. It was greater, however, than the value peromes using the 
 Felation \ given Sturm and Moore, (9,10) 


Pa“ 


for 
; for the of load to8 the 
“deflection is 0. "0270", or 0. 007743 Pa 2/D. Formula (23) is given for a plate 
with corners held and edges free to rotate. ag’ this test, the corners were not® 
“held, so that actual deflections somewhat than by this formula 
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i Fringe patterns for rectangular plate, simpy 


‘supported at ends and with free edges. 
fringe patterns are ‘shown in Fig. 29, the ‘moments and deflections 
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se agreement for the deflections, an 
4 q thor (23) agrees well with the exper imental value. 
moment are about 15 percent lower than the difference 
—_ The Moire values for momen are mut 
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ig. 25. 


vod 
Fg patterns for simply square s steel 
plate) with concentrated center load. 


the test procedure followed in several possible 


of error can be listed. Among these are: 
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ply square plastic. 


Lack of precision in loads, 


Variation in plate thickness, ~ 
‘Support conditions other than 
Changes in dimension of photographic film during BIAS and di ne! 
tortion in the enlargements, | 


ie "Relative motion of model, “grid, or camera between exposure 


using a relatively simple test setup, have successfully kept most of these at Jj 
alow level; however, certain suggestions can be made for further improve- | 


. The greatest source | of error is believed to be in the calibration of the 


naa view of the results obtained, it appears that the procedures followed, 


_ material. For any extensive test program, the use of metal models should be 
seriously considered. Such use would probably require a more rigid support- 


ing frame than used here, ‘since in order to minimize membrane stresses ~ “Aq 


and still obtain reasonable fringe patterns, thicker models with correspond- ; 
ingly larger loads will ber required. The use of silvered plate glass might — 

_ Glass photographic plates may offer a slight increase in accuracy, although 
it seems that the / method used in locating tt the a — ke any error re- 

a sulting from distortion of the film at a low level. 
_ _In reducing the data, the accuracy of the curvature values | epends on the 

- Sn in locating the slope fringes, in plotting the slope | curves, and in| 
determining the curvatures. It seems possible that two different workers, 

F starting from the same set of slope pictures, might obtain moment results == 
varying by | as. much as ten percent, the deflection results , however, vine ob- 
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The results obtained information can be q 
ee. 84 a _ Obtained by the use of Moiré method with the setup described and using plasti 


FLAT PLATES 


Fig. Moments and eflections for square, simply supported 


Pilate w with ith concentrated center load. 
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‘models; with the suggested revisions, precision should be 

3 The method should be particularly applicable in the study of unusual shaped 
bs plates with boundary conditions which are difficult to formulate ae tt 
_ In some cases, however, even though formal mathematical solutions are 
_ not possible or practical, numerical approximations, such as those ad 


with finite differences, may be quicker and more 


"intervals could be sought by superimposing the unl negative on 
“Moments a: and deflections in edge supported ‘plates. ‘combination of the 


- Moiré method and differences could be used in 1 studying the behavior of 
plates with edge beams of various stiffnesses. _ 


& Plates with large deflections. Plates with deflections of varying magni- 


— * _ tudes up to multiples of the thickness: could | be studied for various edge 


ae ‘?~) Plates with reinforced cutouts — a study of the effect on deflections and 


Studies of plate elements in composite structures. 

These are but a few of the | possibilities for further work. In its | present q 
r the Moiré method should yield useful results in these studies, and with 
the improvements in the technique which should come with further and wider 
use, the method, applied separately and in conjunction with other analytical _ 
and experimental valuable information for appli- 
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é ‘spectrum curves which show the effect of s size and ‘distance of earthquake sand i 
: period and damping of structure. High elastic stresses are developed during 
be 


ea ground motion, so the desirability of a type of limit design is indicated. 7 HY 


response of a is affected by the size and distance of the 
and by the natural period and damping of the structure. It is found that during 
: _ strong ground motion ordinary structures will be subjected to o larger geostera- 
§ tions than are prescribed by building codes and, therefore, the behavior of the — 
_ structure depends strongly upon its ability to dissipate vibrational energy 
when it is overstressed. This indicates the desirability of a type of limit is 


_ Note: Discussion open until March 1, 1960. Separate discussions should d be submitted - 
a for the individual papers in this symposium. To extend the closing date one month, 
.. written request must be filed with the Executive Secretary, ASCE. Paper 2220 is 
pat of the copyrighted Journal of the Engineering Mechanics Division, Proceedings a 
ta of the American Society of Civil Engineers, Vol. 85, No. EM 4, October, 1959. - ¥ 
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damage in Mexico city (1) by the earthquake of July, 1957 was 
sharp reminder of the importance of aseismic design of structures. The con-§ 


clusions drawn from observations in Mexico City agreed, in general, with thos lw 
7 derived in United States cities after earthquakes and pointed up once again the | ie 
fact that conservatively designed structures survive with at most minor dam- “a é 
7 age, whereas structures not designed to resist earthquakes may collapse, as a 
; Pan a | number of them did in Mexico City. ‘Structures unconservatively designed £ 5° 
against earthquakes in many instances suffered extensive damage such as 
i= _ Shattering of walls, buckling or fracturing of columns, and failing of beams in sia 
bending. In some cases large multi-story buildings were so badly damaged — e 
it seemed a miracle that they did not collapse. Although the intensity of 
_ ground motion in Mexico City was not high by California standards, it was 3 ie 
appreciably more intense than most Mexico ‘City engineers had anticipated. 
y _ From the point of veiw of what can happen when a city is shaken more severe} Wrsie 
the building code de anticipates, the Mexico City experience was particularly 
7 It was already known from United States experience that improper design 
ae and construction leads to damage during an earthquake so that this aspect of 


the Mexico City experience was not particularly informative. . However, very 
_ interesting building behavior in Mexico City was shown by those buildings that 
i[™ _ were damaged more or less seriously but did not collapse. At first glance, — 


many of these cases seemed paradoxical. For example, some structures wer dista 
= so badly damaged with failed beams, columns, and walls that from a stress F trate 
J analysis point of view it appeared that hardly any lateral strength remained. F 
Since the duration of ground motion was close to 60 seconds, before the end of th 
- the earthquake these structures must have been in an advanced stage of dam- & moti 
” In other cases a damaged building remained standing even though a standarg§ guri, 
lateral strength analysis of the building in its damaged state would indicate > By 
# that its ‘strength | was so ‘small | that it ‘should | have collapsed. On the other ha a surv 


- was particularly striking to observe that the old masonry structures ins 
-— cases came through with little if any damage, whereas the relatively 
_ modern multi-story buildings in many cases suffered severe damage; and this 
was not because the old buildings were particularly v well designed or 
Een Many | of the modern buildings were seriously damaged, whereas 

some om with, low lateral al strengths su with no 


a 


resist a specified static lateral ieee. To understand such behavior, it is — 


necessary to examine the performance of structures during earthquakes from 


"motions. If sufficient were known about the physical properties of the struc- 
tures and about the ground motion during the earthquake, a thorough vibration 
analysis would, of course, completely explain the observed behavior. Al- = 
_ though such detailed information is not available it is possible, by using th the 


results of analyses of past earthquakes to explain the observed behavior. a 


— — 


Earthquake Ground Motions. A typical strong ground motion record is 
hown in Fig. 1. This motion was recorded approximately 30 miles from the 
enter of the fault. Ground motion more intense than this has been recorded 


ring other the > most being approximately 
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sign FIG.1- ACCELEROGRAM FOR. TAFT, CALIFORNIA 
ctof § EARTHQUAKE OF JULY 21, 
ice, . In Fig. 10 is shown the ground motion recorded during this earthquake at a 
S Wel distance of approximately 70 miles from the center of the fault . This illus- 7 5 
ess trates how the intensity of motion diminishes with distance. 
ined. Fig. 2 shows the recorded ground motion at a point close to the epicenter 
end of : of the small earthquake at San Francisco, 22 March 1957. ‘The: duration of the fee 
dam- @ motion was very short compared with that of a large earthquake. This motion 
_ was recorded on rock whereas the motion Fig. 11 was 
tandarg} during the same earthquake on 120 feet of alluvium. 
cate | _ The foregoing ground motions were recorded on U. Coast and Geodetic. 
ha Survey strong-motion accelerometers having a natural frequency of vibration 
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2. Vibration of ‘Simple Ifa a ‘simple s structure, consisting of a 
m with spring constant k and damping c as shown in Fig. 3, is subjected to” 
earthquake ground motion it will be excited into vibration. The structure wi 
- start from rest and will build up oscillations which then will die away at the 
= the mass, and damping of the structure, and will, of. course, 
7 ‘different for different ground motions. To explore all of the various possibil 
Zz s of structure and ground m motion would be an extremely big job; however, 
engineering purposes the interest lies not in the precise form of the vibratic 
but rather in the maximum displacement and stress experienced by the stru 
ture during the earthquake. This information can be exhibited in a conveniet 
form by the following procedure. If sucha quantity as maximum displaceme 
is computed for all different combinations of mass, stiffness, and damping, | 
«this information could be plotted as a type of influence-diagram for the grou 
motion and the maximum response for any structure of this be re 
the diagram. This has been done as described below. 
Newton’ of motion for the systemis 
The e relative ‘motion y is thus the same as fora structure resting on an 
+ immovable base and subjected to a horizontal for force (- -mzZ). As is well- knows 


pen, 
q 
~ 
| 
| 
the 
—— 
EGREE OF FREEDOM STRUCTURED 


vever, 
ribrationg 
ie struc 
nvenier 
laceme 
nping, | 
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ay wadamped of vibration tas 


= time at which y is evaluated 


behavior of the system is made ‘clear introducing the 1e following 


oh. 


Using it can be that the motion of the system is 


described by the following equations for relative displacement and Bin ae 


otential ener of the system 


The Kinetic energy is defined to be 4 m and it is is given by 


total energy is then 


+ KE) = mE(1 + sing. wt - a ) 


The ground motions s recorded during strong earthquakes es exhibit the pi proper- ia 
ai of random functions?» 3 for which E, a, 8, y; are slowly varying functions — _ 
of time, and the motions of: y; ¥ closely approximate a sinusoidal shape. For 
_ moderate amounts of damping the factor (a - 8) may be set equal to zero and a 

‘the he maximum total energy is om by 
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Where is the maximum total energy per unit mass the structure 


during: the ¢ earthquake. ‘The maximum relative displacement, relative velocity 
and absolute acceleration attained by ‘system may be 


a Vv 
the notation ‘Sy ist is used for 3 The” maximum Ss) shear for force 


The maximum inertia force acting on she: mass is 
which, of course, agrees with an). wa. 
_ The maximum response of the structure is thus characterized by the "7 
Sy which for moderate damping may be 


gil It is — ‘that the response of an n undamped structure builds up roughly | 


proportional to the square root of the duration of the earthquake ground 
x motion, whereas, the response of a damped structure builds up until the ener- 
gy loss balances the energy input. 
_ A one-mass system that can vibrate | in both nestennted directions will re- 
spond to both horizontal components of ground motion and the resultant motior 
_is given by the superposition of the two component responses, but the values 
of S,, will not, in general, occur at the same time for the two components. 
3. “Complex Structures. A multi-story building will have several modes of 
_ vibration and consequently its response | to ground motion will be more com- Bf. .. 
J plicated than is the case for 2 a simple one-mass system. The general nature § 
AN ood of the motion can be elucidated by the following considerations. Consider a 
general elastic structure vibrating in the X-Z plane as shown in ‘Fig. 4, - Let . 
4 the structure be linearly elastic and let it be. such that its free vibrations cang —s 
resolved into normal modes of vibration that are orthogonal. 


4q 


where 


incr 


4 = = shape ofi- th node of 


‘This follows from the fact that earthquake » ground acceleration is essentially ae 
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FIG. 4- ~ MULTI-DEGREE OF FREEDOM St RUCTURE 


As in the case of nee oscillator, the response y y to base acceleration — 
ef a is the same as the response to a lateral force -mZ acting on each element of 


mass sm. The effect of Zz acting for an increment time dt is to produce an 


and from the orthogonality of the oad 
dm 
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where dm that the are to be | evaluated over 
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om the case of the se oscillator, the : response mse to a continuous Z is s the 


™ 


‘The expression in the ee is ak aaa same as s equation (3) which gives the 
% response of a simple structure having a natural undamped period T oi” a 
response of the complex structure can thus be determined if the mode shapes 
_ are known and if the responses of simple structures of the appropriate periods 
anddamping canbe computed. complex 
computation of the complete response of a complex : mructare is 
extremely laborious and is best done by computing machines. — 


il 
Similar expressions involving Sy can be written for the maximum forces, 
_ stresses, shears, bending moments, etc., if in each case the maximum — tides 
absolute value for each mode is added. (4) The analysis of structures more 
complicated than the above can also be carried out in a similar fashion and % 
results similar to equations (14) and (15) will be obtained. 
e _ From equation (14) it is seen that the true maximum response will, in od 
general, be somewhat less than the ‘Tight side of equation (15) since each term 
positive and negative values, and in addition the amplitude functions, V2Ei, a 
vary with time. Thus, there is only a small probability that all of the terms 
will reach their peak values and be additive at exactly the same time. Studies 


ground motions, if only two well separated modes of vibration are involved, 


equation (15) overestimates the true maximum by about 5%-10%. However, if § 


three or four modes are involved and the maximum of each has approximately — 
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ie sec ‘made at the California Institute of Technology show that for strong earthquake © a 


- MAXIMUM VELOCITY - FEET PER SECOND 


the same value, equation (15) may significantly overestimate the true maximum q 


value. - For such cases it has been proposed(3) to use the geometric mean of the 


4. Earthquake Spectra. ati 
— S, will describe the maximum response of a wide class of structures. The 
is called the spectrum|(5, 6) of of the ground motion and a pot of 


shown by equations (10), (11), (14), (15), the values 


Sp, whic 
_ _ 
— i “It th is damping in the structure e 2 ode | 
| 
|, = 


EARTHQUAKES 
in effect an influence diagram for the maximum a simple structure, 


Sp which is equal to Sy; is called the displacement spectrum and 


picture of the 1e effect of of the +o motion and in ‘this anes it characterizes the "4 
+ A plot of Sy for a typical strong earthquake ground motion is shown in Fig. 
5. ‘This gives the maximum velocity response over a significant range of por 
periods and for 0%, 10% and 20% of critical damping. The shapes of these 
curves are typical for ground motion measured moderately close to the ae 
jepicenter of a large earthquake. In fact, when the normalized spectra for 
jeight components of the four strongest ground motions yet recorded ‘€l Centro 7 
1934, El Centro 1940, Olympia 1949, Tehachapi 1952) are averaged, a quite A 
smooth set of curves is obtained. The average velocity spectrum for strong e : 
ground motion shown in Fig. 7 and the corresponding average acceleration | Sie 
Sspectrum shown in Fig. 8 are based on the average of the four shocks. eo 
_ The velocity spectrum for the ground motion recorded in Hollywood during © 
the 1952 Tehachapi shock(7) is shown in Fig. 6 and this is very similar Rout * 
that for the ground motion recorded at Pasadena, 4 both | recordings being a j 
approximately 70 miles from the center of the fault. It is seen that the shape e 4 &§ 
of this spectrum is different from that of Fig. 5 which is for the same earth- a> 
iquake but at approximately 30 miles from the ‘center. The short p period end of - Wi 
the spectrum in Fig. 6 is relatively depressed as compared to Fig. 5 thus Ci a - 
indicating that the higher frequency components of the ground motion attenuate : 
with distance more rapidly than the lower frequency c components. ae oe 
_ For ground motion of relatively small earthquakes, such as the shock at 
San Francisco, 22 March 1957, recorded quite close to the he epicenter the 
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E-W “VELOCITY SPECTRA, , HOLLYWOOD 


"velocity spectrum is found(8) to be relatively high in the low period end. . Fig. E 

* 9 illustrates the differences in shapes of spectra according to whether the JJ 
shock is large or small and whether the epicenter is near or distant. . The ay 
curves represent approximately the Taft, 21 July 1952, Hollywood, 21 July © : 
1952, and San Francisco, 22 March 1957 spectra. igge is 
‘The | ordinates of Figs. 7 and 8 ‘Should be multiplied by the 


El Centro, 18 May 1! 194¢ 


Taft, 21 July 1952 
. Vernon, 10 March and q 
Spectrum Intensities. The maximum of structures are indicat- 
ed by the S, spectrum curves and, hence, the area under one of the curves is 
7 a measure ‘of the intensity of the ground motion in the sense that it expresses 
a _ the average response of structures of that amount of damping over the range 
periods. The spectrum intensity (9) is defined be the 
curve between 0.1 and 2.5 seconds period: 


In this expression n indicates the amount of damping. It Aina te nebes that 
the spectrum intensity is an objective measure of the maximum stresses that 
Hay Ww would be induced in elastic structures by the ground motion. On the other — 
Cm ha and, the comomnly used Modified- -Mercalli intensity is a subjective assess- 
ment of the degree of damage actually by the ground motion. 
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FIG. AVERAGE VELOCITY SPECTRUM CURVES 
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ee for the two horizontal components of motion and also lists the aver- 

age intensity of the two components. The average for the two components of 
. the 20% critically damped spectrum intensity is also listed. . The earthquake 4 
-Magnitudes (M), which are measures of the sizes of the earthquakes, are also 
_— listed. It should be noted that earthquake No. 4 (July 21, 1952) was produced © 

by more or less vertical slipping on the fault rather ‘than the more usual 
horizontal slipping and the motion at Taft was not so intense as might be “4 
expected from a Magnitude 7.7 shock at this distance. Records” No. 9 


foundation conditions and hence are not directly comparable to the other 
ay. _ It will be noted that the ratios between the zero ro damped samen: are 4 
not the same as the ratios between the 20% damped intensities. This is 
_ because the damped intensities are measures of only that portion of the re- 
cord that had the greatest accelerations and the duration of the ground motion 
does not influence this intensity. On the other hand, the zerodamped 
ae intensity is a measure of both the magnitude of the accelerations and the 
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, Measured Building Motions . The foregoing analysis has been checked with 


_ ground motion been recorded together with the resulting building vibration. In 
“these cases, ‘the application of | the ‘spectrum method to computing the maxi mun 


_ Fig. 10 shows the recorded aceleration 10) of the roof of a 10- -story rein- 


; "measured building motions although only in a few cases has moderately strong " : 


IN concrete building and also shows the earthquake ground motion record: 


oe ed in the basement.(10) It will be noted that the maximum ground acceleration}. 
wt 5% and the maximum recorded roof acceleration was 14%g and came most 
ly from the first mode of vibration. Fig. 11 shows the recorded acceleration 
on the 16th floor of a 16-story steel frame office building and also shows the 
ground motion mn recorded in the | basement. . ze It will be n noted that the maximum 
= ground acceleration was 5% and the maximum recorded 16th floor accelera- 
_ tion was approximately 12%g. For this case, analysis showed that the second 
_ and third modes of vibration contributed strongly to the building acceleration. 


Fig. 12 shows base and building accelerations recorded during an n explosive- 
generated ground shock. 
It is seen from the that appreciable building accelerations are 
produced by only moderately strong ground motion. In the case of strong 
ground motion such as the El Centro, 18 May 1940 record accelerations of the 
order of 30% to 40% are indicated for buildings 1 like the Hollywood structure. 
; Most buildings could not withstand these accelerations without being o over- pe 
stressed. On the other such structures have survived rong ground 


Vibrational Energy of Structures. ~The undamped, multi - -degree of freedom 
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FIG. BASEMENT AND ROOF ACCELERATIONS, HOLLYWOOD 
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and equation (18) 1 may thus be written(12) 
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Although: the average zero damped spectrum curve is not a tiesiieaaiad line 

|e is seen from Fig. 7 that in the range T>0.3 seconds it is not very far from - 

it t. Accordingly f for structures ‘whose ‘pertinent modes have above 0. 


depends only upon the total mass and on and does not depend upon 
how the mass is distributed or how the framing is laid out. 1 afd on on : 
ry . Effect of Damping. It is seen from Fig. 7 that structures with low damping 


states that on the average the maximum total energy fed into the structure | = 


:, will be subjected to relatively high stresses in the event of strong ground = 

motion. “Moderate | amounts of damping will give appreciable reduction of 

— stresses. It is also clear that large amounts of damping will give large re- ‘ i 

f 


_ductions in the forces and stresses in a structure. Vibration tests show that _ 
for small vibrations of ordinary structures the amount of damping is of the SS. 

order of 10%of critical which is not large. Structures of the types of the __ 

Hollywood Storage | and the Alexander buildings could experience accelerations 
= of the order of 30 to 40% during : strong ig ground ‘motion if they had 10%damp- | 
ing. However, if such a structure is subjected to strong vibrations, the 
consequent cracking, rubbing, plastic yielding, etc., , can absorb large amounts — 
of energy and effectively give the structure large damping. Most structures _ 

have considerable ca acity for energy dissipation when they are overstressed. ' 
An exact analysis(13, 14) of the vibrations of overstressed (non-linear) ae 


structures is ; exceedingly difficult and the results of such analyses ; are ‘not = 


Poa 


yet in a form to be used for design. An approximate method, — on er: 


tion (19) can be utilized. ' (1 5) 
Cary) 


Consider a structure starting from rest and excited into vibration strong 
oot motion. The vibrations of the structure will continue to increase in ra q 
amplitude and eventually some part of the structure will be overstressed. ‘fl —— 
the ground motion is sufficiently strong the amplitude of ‘vibration will con- 
tinue to increase and overstressing will occur on each swing of the structure. 
‘After a time, the intensity of the ground motion will decrease and the ire 
will quiet down to the point where it is no longer overstressed during the a 
and following this the vibrations will | be gradually damped out by 


the ordinary damping inherent in the structure. This process can be describ- 3 


—_ E represents the total energy input into the i la to the time that 

the last non- -linear overstress occurs; D is the energy dissipated by the _ ae 

structure up to this time: is the potential energy inthe structure at the 

time that the last energy- absorbing overstress occurs. 
If, as in equation (19), it is taken that 


he 


-£ this assumes that the input energy into the overstressed structure is tl the tae 


| game as the maximum input e energy for an equivalent, elastic, ‘undamped 
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October, 1959 
‘This equation specifies: the amount of energy that must be dissipated if the 
ashen is to remain standing . For ground motion of approximately the — 
intensity of the El Centro 1940 records the equation may be written 


_ perties of the structures and the intensity of the ground motion to make an © lative 
analysis, equation n (22) agrees ‘reasonably well with field observations of with 
a The e importance of plastic energy absorption in controlling vibrations is 
illustrated by the behavior of a steel bar of 1 sq. in. cross-sectional area and 
Pie length. If the elastic limit of the steel is 33,000 psi the strain ae 
= (0. 0.0011) = 1.5 ft lb per cu in. 
_ If the bar is stretched plastically to a strain of 0.02 the ‘energy absorbed © 
plastically is (assuming perfect plasticity) 


33,000 (0. - 0, 9.0011) | 50 ft lb cu ‘ 


ings 


i the yield point and could absorb an appreciable amount of Ihsosth and still | 
have an adequate factor of safety against collapse. For example,(15) a rod- 


_ acceleration to remain elastic during ground motion of the El Centro 1940 
intensity could be designed for 12%g and undergo plastic , deformation and 
- still have a factor of safety of approximately 2. A plastic design must, of 
_ course, be done very carefully to ensure that rupture or brittle fracture does 


iz 7 braced, elevated water tank that would have to be designed for 40%g lateral 


+ 


‘not occur. The total energy absorbed in failing members by oscillatory 
_ plastic strain is not as yet well established so that the the Gosign should be 


— those apparent anomalies that are frequently observed following a> 


ae a) After almost every strong earthquake it is observed that a few old ll 
ian. apparently very weak buildings survived even though their static later — 
- strength (as computed on the basis of a a Fg load on the structural frame) was § 
low. This can be explained as the result of two factors First, the 
A true ultimate lateral strength of a structure is often much larger than the — 
‘nominal design lateral strength. Sec ‘Second, such old, loose-jointed structures © 
may: absorb large’ amounts of energy during only moderate vibrations and, gq 
= is seen from the spectrum curves, such a structure could survive strong 
motion without actually being subjected to high stresses. 
It is not uncommon to observe that the same (ground motion that 
ee apparently did no damage to an old weak building damaged a relatively stroll 
building. This is ‘readily explainable. Consider a well- designed and construc: 
oe ed building whose lateral strength was so high that it could resist a static Hy 
en load of 25%g without being overstressed. If the design and constructid 
_ were ‘such that the structure had only small damping during elastic es 
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it wa 
— veight of the structure. For such ground motion approx-§ | 
a a 4 mately 1/10 of a foot pound of energy must be absorbed for each pound = the 1 
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some member would then see. distress. If a designer were sufficiently 
clever at incorporating damping in structures he could design a 1 weak struc- 


ture that would be undamaged and a strong structure that would be damaged ‘#4 id 


vy ¢) ‘After | the Mexico City earthquake, whose epicenter was 220 miles away, 
it was observed that in the same district old masonry buildings came through — 

- I the earthquake on the whole better than modern steel and concrete buildings — aa 
ppror- | that had greater lateral strength. This can be explained as resulting from — e a 
Of Fi the long 1g distance of the epicenter from the city. As illustrated in Fig. 9, ata a 

~ § long distance from the epicenter the short period end of the spectrum is re- 
latively depressed. old masonry buildings were mostly low structures 
with largely multi- -story, long-period buildings, so that in effect the old build- _ 
ings experienced a much weaker ground motion than did the modern buildings. 
If, in the future, the epicenter of a relatively small earthquake, : such as shown 
in Fig. 2, were to be located close to the city, it could be expected that the © 

} low and stiff buildings would suffer more than the tall , flexible structures. we 

After almost every destructive earthquake some buildings are observed 


to be in 1 such a badly damaged condition that, from the static design point of 7 
view, it seems a miracle that they did not collapse. It is often said that if 
the ground motion had been but a shade more intense many buildings would ; 
§ have collapsed. If the damage to a structure is looked at as a process for 
absorbing energy, the survival of these structures is not so mysterious. In 
these cases it is not so much a matter of static lateral strength as it is a 4 . 
question of whether | the structure c can continue to absorb sppreciable a 
» during its vibrations. . For example, consider a badly damaged building that 
_ @ has cracked a large nu number of interior partition walls, has extensive crack- 
ing in exterior masonry walls, and exhibits many cases of severe overstress_ 
inbeams and columns. The very fact that the damage is so extensive shows | 


that a large amount of energy has been dissipated and it is likely that further  : 


vibrations will continue to dissipate energy. ‘The structure may thus not be — 

in so perilous a condition as would be concluded on the | basis of a static ere “a - 
al force analysis. . Opposed to this, consider some one- -story concrete- -shell 4 
structures that collapsed in Mexico City. These were long, narrow concrete . 

| roofs, each supported on two concrete columns. . With a heavy roof supported ‘ey 
on only two columns the inherent damping for small vibrations | would be seat 
small. Furthermore, } plastic energy absorption could take place only in a 
|overstressing the two columns in bending and this was not very large. ‘iso, 

one of columns lost its its strength the entire ‘Structure must pecegsarily — 


really more surprising than eT some of the — 

10. Conclusions. ‘The response spectrum cu curves, which are a type of influ- 

ence diagram for earthquake loads, describe the important characteristics of _ 
strong earthquake ground motions by exhibiting their effects on —~ 
The spectra also provide a logical measure of the intensity of earthquake r 
ground motion which has engineering sugnificance. By means of the —— 


curves the apparently ; anomalous behavior of some structures during ¢ — 


construd _ If the point of view is taken that structures should be designed to remain PP * 
atic plastic (without damage) when subjected to earthquake ground motion, Sook .. 
structio satisfactory procedure can be developed using the spectrum curves. s. It is 
bratiol that be more than the used equivalent 
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static but it would be so. real 1 difficulty a arises 


ui region all buildings should be designed to survive the strongest possible a 
ground motion without any damage. it is more reasonable to take the point of 

. view that the design should be such that the structures will survive the more 
7 frequent, moderate ground motions without damage but in the more rare event 
of very strong ground motion, damage would be tolerated so long as it was 

' not a hazard to life and limb. This is the usual point of view of engineers in 
earthquake areas. It is seen from the material presented in this paper that 
,_ this point of view is equivalent to requiring two ‘different design analyses for 
= “ each structure. One would be to ensure that no ‘damage resulted from mod- 
f 4 - ‘a erate, elastic vibrations, , and the other would be a “limit design” type of 
a analysis to ensure that there was sufficient energy- ‘absorbing capacity to © 


F te give an adequate factor of safety against collapse in the event of extremely 


it is s not difficult to ‘make | a 4 reasonable limit ines along the lines set forth 

in this paper. However, in the cane of more complicated structures the | 
_ problem becomes | very difficult, particularly because we do not at present 


_ of ordinary building materials and structures during vibrations. This is an 
_ important problem now facing engineers, and it is being attacked k by research 
~ projects sponsored by the California State Division of Architecture. One of bi 
these projects is a study of how to compute the ultimate lateral strength of 
ee 7, structures; another, being of structures; a third, being | done at the California 
qj Institute | of Technology is the development of a large- capacity shaking ma- - 

_ Chine with which it is planned to vibrate some typical structures to the pois 
of serious damage and measure their energy absorbing capacities. 


projects | are re just a beginning, and more research is needed on the. general 


tt is, of f course, true that the lateral forces prescribed fc for r design by build- 
ing codes. are based to a certain extent on observed damage behavior of 


_ structures during earthquakes and in this eense they do reflect the energy 4 


requirements to non- typical and unusual structures is not warranted for the: 


energy absorbing capacity of such structures may be quite different from 


id absorbing properties of f typical structures. application of these code 2 
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“PLASTIC 


In the application of plastic theory to design, just as in the case of allowa- i 
ble stress design, there are a number of factors the importance of which n must — 
be evaluated. Those considered in this chapter include shear force, local | —_—_ 
buckling of flanges and webs, lateral instability, and repeated loading. The 
problem of of column buckling is treated il in chapter 7 7 (Progress Report _ 


m 


_ Simple plastic theory is based on the concept of the attainment of the full 


plastic moment at certain sections, followed by adequate rotation at this 
constant moment value as the applied load increases. It must be recognized, 
however, that the magnitude of the moment at which large inelastic rotations 7 
occur is not constant under all circumstances. In calculating the full plastic ie 


moment Mp as Z, it is assumed that the member is stable and is eae = 


to pure bending; thus axial thrust and shear force are both ignored. 4 ve Ae 

+ Since plastic hinges usually occur at positions where shear and axial force 
are present, it is of importance to be able to predict the changes in the values Y 
_§ of the plastic moments due to these causes. . In most practical cases the in- 
fluence of shear force will be very small. Ps In some wpectat cases the combined — 


 Horne(6-1) considered the effect of transverse shear force for the case of 

a beam of rectangular cross section bent about one of its axes of symmetry. 

An approximate solution was obtained, valid for shear forces less than a cer- - 
tain limit. This solution was also extended to an I-section bent about the Pa. fa 
Leth(6. 2) modified this latter solution to dete 
value of the full plastic moment in the presence of shear force, valid over the 
full range of shear forces. He included an upper bound on the full plastic mo- 
ment when the shear force has its full plastic value. Heyman and Dutton(6. 3) 
proposed a semi-empirical theory which is also valid over the full range of oa 
Shear forces. -Green(6.4) derived an upper bound on the value of the full plas- 
tic moment in the presence of shear force for beams of both rectangular and — $ 4 
I-section. This upper bound agrees closely with Leth’s lower bound for io 
g Values of the shear force. Upper bound solutions were also derived independ-— 

ently by Green and by Onat and Shield(6.5) for the case of a beam of rectangu-. 

lar cross section, treated as a problem in plane strain : and assuming tl the mate- 7 
tial to be plastic-rigid. Drucker (6.6) demonstrated the lack of uniqueness of 
interaction curves for shear and moment, explored bounds for a beam of | rec- 
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* sania 7) derived an expression valid over ar te e full ra range nge of of shear, in gen- a 
eral terms, which is applicable to both rectangular and I-beams. _ Fujita(6. he 
' - developed a lower bound solution of a cantilever beam with rectangular cross 


a section, assuming the plane stress state and allowing warping at the fixed end. 
- This solution was modified for wide flange sections by Beedle, -Thurlimann 4 


S -_ Ketter(6- 9) including an extreme case of shear failure. be: The effect of both B 


results obtained in Ref. 11: are in fair agreement with 


Analysis 


In most practical situations, or sections are used. the es- 


sentials of the problem can be ‘shown by the simple « example of a cantilever 
beam with rectangular cross section of depth d and width b (Fig. 6. -, The 
following assumptions and limitations are made for the sans acne 

1. The problem : is analyzed under the assumptions made | in pie trea aa 

theory as a problem of plane stress. 

2. Only the on portion can carry shearing forces. 


+ 3. The Mises - Hencky yield. criterion is adopted for the combined cred 
4 


aa; = - The fixed end is allowed to warp. 


-s the beam and start to penetrate into the cross section — 6. my _ External 
and internal moment and shear at section S-s -s in Saal (6. -la are given b “a 


y 4 


_ih order to determine the shearing stresses > the horizontal equilibrium 
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1. wee ‘Zones d Distribution is a 

any shearing s stresses T, since the resultant 
forces on both faces are equal. Hence ‘no shearing ‘stresses exist in the plas- 

tic zone. The shearing stresses” % in the — axis along section dl Bre- 
sist the unbalanced horizontal ad 
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4 
tt is assumed that ies resistance of the section will be citieitnd if this shear- 


stress reaches its yield value at this centerline the 


% wih “yield 
‘The limitin limiting value value for is given by Eq. 6.7 as 


nies oe condition will first be reached at the fixed end x = = 0; therefore 


Replacing y Eq. 6.9 
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PLASTIC 
where modulus of rectangular cross section. 


6. 11 the desired relationship between the plastic n moment modi- 


fied by y shear and the ene plastic ic moment Mp. It can can readily be § solved for the oi 
parameter a/d (le ngth of cantilever/depth of section) for different values of _ 
the ratio ps/Mp. _ The result is owe in Fig. 6.2 for the curve marked 


ate plastic eatibes Z of the wite- flange section is used when yielding due to 
normal stresses has penetrated into the web of the section and when the web— 
thickness w is used instead ofthe breadth, b. 
16 Fig. 6. 2 shows the comparable reduction of the ful full plastic moment due aan 
shear force for various ratios of the flange to web areas Ag/Aw- _ The case 
At/Ay = 0 ‘corresponds to a rectangle. The | dimensions are given in Fig. 6. . 
As indicated on the inset sketch of Fig. 6. 2, only that portion of the chart for 


” ih 0. -65 is st shown | to the large scale. 
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Note : | 
which does 
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web in shear. 
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Fig. 6 2. F Reduction of Plastic Moment by shear | Force ies > 


6.8) 
6-75 
oll | 
shear- 

(6.8) 
— 
(6-9) 
efore — 

hie 


Itis ‘emphasized that this is a maximum ca carrying but 
- eather corresponds to “first yield in shear.” It would therefore be expected — 
that members tested nogeedl this an would display a greater | ‘strength ‘tha 


cases of high te values an extreme condition may de- 


velop where | the entire moment is taken by the flanges | and shear produces | 
yielding | oft the as indicated in Mig. 


: 
Using these two relations the 


(6.12) 


‘The curve for complete of the in (Eq. 6.12) shown 
in 6.2 was for an average ratio 


An approximation to condition full Ul yielding of the web 


there obtained for the maximum allowable shear at load, 
Vy" = 18,000 wa 
where ‘Vis in pounds and w and are in inches. 
2 ‘maximum shear force at ultimate load V is greater than 18,000 wd ow 


a cas the web could be reinforced or Eq. 6.12 could be used to predict the modified 
plastic moment value. When is less than V Eq. 6.11 is appropriate; 


however it will be shown later that the in moment 
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6.3. Assumed Stress Distribution ina WF 


d 
‘Section Loaded by Shear and Moment 


7 i, in addition to shear, an axial force is also present a similar ee of _ 
t analysis may be applied when both : shear and axial force are present. For | 
wide-flange sections with At/Aw = 2.0 and dg/dy = 1.05, the problem has has” 


i } been solved i in Ref. 6.11. Fig. 6. 4 shows: the relationship of Mp m/Mp to to = | 


a/dw for a cantilever beam, the parameter being the a of axial load aa 

plastic hinge moment n modified to inchete of effect of of 

= compression and shear force 


- = 0, is the s same as’ that given in 


7 carve for axial force equal to zero 


> 
2 = 2.0. other curves sho the nce of 


‘occur in combination at points where plastic hinges are expected. . For such — 


an eventuality, however, Fig. 6.4 may be used to modify a trial “design that oe 


_ Tt is not easy to make a precise comparison between the theoretical pre- a 
dictions and experiments, since the effect of shear is quite small for wide- 
flange sections until the shear force V approaches the extreme value of © 
| Shear failure Vy as given | by Eq. 6. 14. When V is in the neighborhood a I 
y the load- cu curves for beams show no ) appreciable horizontal por- 


tion corresponding to plastic hinge action. - Instead, it is found that the load-— a | 


pounds), deflection « or the moment-rotation curves continue to rise fairly steadily and 
codified Such an extent that the full plastic moment is exceeded; henceno 


capacity m However, the intersection of two tangents to the load-deflection or moment- 


the and range provides one possible criterion 


— 
4 
(6.12 
2) sh 
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f the web ~ 
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‘Fig. 6.4. Moment eatin as Modified Combined 


‘brik 
Action of Shear and Thrust on WF Section 


for evaluating experimentally the modified plastic moment of the section due 


to shear force. Beyond this value of moment the deflection or curvature of 


the beam increases more rapidly rer eet 

jan The Gctermination of aps is ‘shown in Fig. 6.5; the t test data is taken from 
‘Ref. 6. - Predictions from Eqs. 6.11 and 6.14 are also indicated for | com- 
parison. The shear force V in this: case is greater than 18 ksi, Eq. 


4a ‘Fig. 6.6 shows in ‘tabular form a number of tests which indicate the influ- 


at left. is indicated the size of member and the 
; p the depth of the member and “ “a ” is the ratio = ion constant shear, “a” is” 
the ¢ distance from a} a plastic nina to the point a inflection. - Thus, the quantity 


“a” is sometimes referred to as the length of an “equivalent cantilever beam: 


4 To the | right is | a bar graph on which is plotted the maximum moment corres- 


due to Shear, the latter being Gcermined The 


| 
~‘Mmponding to the observed Maximum load and observed reduced pla 
iii q 


ction due 
ature of 

-aken from 
or com- 
ce 

the influ- 
shown 
is 


“a” 
quantity 


ver beam? 
nt corres- 
c moment 


é . theoretical 1 limit according to Eq. 6.11 or Eq. 6. 12, (depending on which came 


web Vy The exception is Beam | No. | 13 which was 0 deficient 


f utilized because of practical limitations. A beam with pure m moment, for ex- 


} cannot be counted upon to support a moment greater than Mp 


- 


‘Fig. 6.5. Determination the as Modified by 
Shear, Using Tangent Intersection Method 


is the controlling one) is shown by the dotted line. In the range where the 
shear force V is high and greater than 18,000 wd (pounds), Eq. 6.12 is used 
to compute the limit predicted by the theory. ‘This figure shows that the 
actual strength of the structure usually exceeds the full plastic moment Mp 7 
if the shear force Vv is less than 18,000 wd (the value which assumes full © . 


“exists. Even when the shear force is greater than 18, 000 wd, the test results — 
| show a moment strength greater than predicted by Eq. 6. 12. It is evident, — Pa 


however, that the latter is a reasonable limit. 


ie: While it might ap) appear ar from these results that a more than ample load fac- 
tor exists with regard to shear force when V < 18,000 wd, this fact cannot be E 


ample, usually buckles when strain- -hardening commences, and therefore 


nt The test results on the effect of both shear and axial on 10) s show : a 
ae similar to that of the | bar graph in Fig. 6. 6. e The observed values @ 
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The. effect of shear on the full plastic moment will be for 

frames. High shear and moment occur in very localized zones so that strain-— i . 


ice will set in quickly and, in most cases, permit _ moment to ‘each a 
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Fig. 6. ‘Summary of Test for Beams 
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"PLASTIC 


is 
required for the effect of sl shear if if at ultimate load v is less than 18, 000 wd = 


Combined Moment and Shear Force 


The design is satisfactory as regards shear force 
ov if its magnitude at ultimate load does not exceed i 
| 18,000 wd (V in pounds and w and d in inches). . No ay 
modification of the plastic moment is required. | 


Buckling 


buckling will not occur r for ‘members proportioned in accordance with 
allowable stress specifications. Plastic design requires that a member have 
capacity to sustain sufficiently large plastic hinge rotations so that such — 

hinges may form at certain sections without the occurrence of inelastic in- a 
stability and consequent loss of moment capacity. Hence plastic design _— 


can be treated separately as slaie buckling problems making use of ee 
boundary conditions and restraints pene from the adjacent plate. 


Previous Research 
_ ‘The pr problem of elastic and inelastic buckling below the yield stress is dis- - 
cussed by Bleich.(6-19) pijlaard(6.20) was the first to apply the deformation © : 
theory (finite stress-strain relations) to the plastic plate buckling problem. 1. 
The theory was further developed by Ilyushin(6-21) and modified by 
Stowell.(6-22) winter(6.23) called attention to the problem of local flange - 
buckling, discussing possible limitations in width-thickness ratios. The flow — 
theory stress-strain relations) w: was applied by Handelman and - 
Prager(6-24) to perfectly flat plates. The latter theory does not compare well z 
with experimental results. Onat and Drucker(6- 25) suggested that the devia- _ 
tions between theory and experimental results were due to the — 


| 


4 


fections of the plate. Bleich(6-19) has proposed an approximate method for — 
the analysis of ineteatie local buckling and has considered various structural 


d 
3 
Ap 


Maaijer and Thiirlimann(6: 27) have presented the results of further study 
of the — of plates ‘compressed beyond the yield point and even into the | 


were also presented in this ‘this reference. = 


The theoreticai analysis of Ref. 6. 27 is based on the following assum 
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The material follows an idealized strain 


2. Yielding occurs ‘discontinuously in slip bands such that the material 


3. The in the elastic range is homogeneous and while 


_ inthe strain-hardening range it is homogeneous and orthotropic. — _ 


a 4. A Uinear strain variation is assumed as shown in Fig. 6. 
5. Yielding of a compressed element starts either from | both oe 


_ taneously or from the of the element. 


6. An incremental stress-strain relationship is used. However considera- 


tion is given to initial imperfections of the plate by introducing effective 


7 Interaction between plate elements is taken into ac account bya oe 


_ Assuming that buckling | occurs without strain ‘reversal, ‘the ‘equilibrium | 
for a plate element of of orthotropic material 1 yields: 


4 


ws 


i 


6. Idealized Stress - Strain Diagram for for 


order to derive a rational plate buckling equation for the inelastic range 


between the proportional limit and the point of strain- -hardening, residual | 

stresses must be considered. _ However, no general solution to the problem 

has been developed to date. Therefore an empirical transition curve has_ | 


Te — reader is referred to Chapter 3, of “Mathematical Theory of Elasticity 
by I. S. for discussion of orthotropic material. 
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usidera- 
effective 


efficient 


6.8. Assumed Linear Strain | Distribution (No Strain 
if 
ive 


> 

H = Dxy+Dyx 

tangent 1 in shear 


The reader is referred tc to Table 6.1 1 and Figs. 6. 9. and i 6. -10 for solutions of —_ 
2sidual § various limiting cases of Eq. 6. 15. References 6.26 and 6.27 may be con- A 


ve has. 4 effective plate moduli Dx, Dy, and Dxy and Gt were in Ref. 
‘96.26 for AT-type steel. Their values are as pei 


Elasticity 


| 


an 


‘Table Solution of 
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EXPRESSIONS FOR CRITICAL STRESSES 
CONDITIONS | AND FOR EFFECTIVE 
: 


lat 
Supported) Oe, = 


y= Fixed 
=b Free” 
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‘the Flange Element of a WF shape 


v4 


= VaDy = 100 


— 
- if 
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Web Plate of a WF Shape 


It should be _ in cies that the foregoing expressions and moduli are appli- 
cable only for A7-type steel plates stressed to the onset of strain- ~hardening. a 
The problem of an axially loaded plate subjected also to bending in the — 
ee of the plate is of practical importance. The web of a WF beam sub- 
jected to both axial load and bending moment presents such a case. This i 


problem has been studied in Ref. 6.27 and a diagram is presented for the : 
allowable — - ratio of the web of a WF section. As shown in Fig. 6.11, the 7 
curves are plotted for different ratios of maximum strain to yield strain, —. 

As the maximum strain of the compression flange increases, the critical ooh 
ratio decreases. The critical value of the axial fo force P also 


"creases as the ratio of the depth to the thickness of the web. increases. 


Wa Tests were carried out to the theory for shapes. 

_ It was found that the theoretical curves gave a good description of the actual i 
of flanges and webs. Fig. 6.12 gives the experimental correlation it 
with theory for the buckling of flange plates. ‘Fig. 6.13 gives the same com- « 

hoon for the buckling of web plates. - The specimens tested in this ee 
included 8WF24, 8WF35, 10WF21, 10WF33, 10WF39 and 
experimental d data are available to with the theory inthe 
j case of bending with a axial force (Fig. 6. 11). Partial correlation | may be — 


_ by examining the results of some corner connection tests, (6.28) if one accepts 
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| | 


“October, 


= Yield ‘strain 

Maximum strain of 4 


compre 


ig} 
“Fig. 6.11. “Ratio io of of Fully 
the fact that the conditions of these tests , were 8 somewhat ‘different from those 


upon which the theory is based. . This comparison is showninnon- 


; ar eae form in Fig. 6.14 where the rotation of a selected length of 


member i is plotted vs. the — - ratio for — = 0.1¢ 10. curve is obtained 


taking cross-curves of Fig. 6. 11. ‘Test not check the 


precisely, but they do indicate a large rotation capacity. = The 14WF30 shape Ff 


an amount less than predicted; on the « other hand, ‘the SOWF108 


greater - ratio a higher load much more 


For flange elements Fig. 6.12 shows that if 17, then 

uniform compression can be strained to the point of strain- hardening before : 
local buckling influences its carrying capacity. 


Fig. that the web of a shape under uniform compression would 


be limited to —_ <33 in order to reach the onset of strain-hardening. ‘How-— 


4 a no such strain is needed in a plastically designed structure and it would 
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x a tests 

Theoretical Curves | 
for o = 36 ksi. 


Restraint 
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‘Fig. 6. 12. Results of Tests on WF Showing Correlation 
with Predicted Strain at Buckling of Flanges 
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Ww Compression tests 


A=Coefficient \ 
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For and compression a design approximation | may 
developed that is based on the given in 11. This 


neon. — 
20 ‘wa — 
>> | — 
\ \ | 
— 
— 
— 
— 
on would on the basis of Fig. 6.13 the limiting value 
g. 
value. 


oO 


. 6. 14. “Results of Tests on WF shapes in which Both Axial - 
Force and Bending were Present. The dotted line 
_ shows the predicted web buckling curve for P/Py 


ine, 


for u use design ‘involving of ASTM -type 


| Ge = 33 ksi), the ) following values are suggested to assure that pre- 


| 

og At = r solutlo Tha 2 

solution sh s: A/Aw the | 

‘a ‘proxi d/w = 
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Buckling of Structural Shapes 


or sections, WF beams of Welded H- of -type 


= 


Uniform co compression 
Bending and 


= 


where b and are width and t of the ulange 
- respectively, d is the sect: section n depth, and w is the web 


oe other outstanding flanges and pl: plate stiffeners in 

tended to be stable in the plastic region. 


For rolled shapes, a small upward variation of the b/t ratio to include, for 
example, the 14WF30 (b/t = 17.58) would be reasonable. This section per- ~- 
formed adequately under test. 
te _In practice, if a selected section cannot meet the above requirements, a 

| Sencees section may be used for stiffeners could be added to prevent local : 
buckling. It is fortunate that nearly all rolled I and WF a sas are satis- 
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ADDITIONAL NOMENCLATURE 
See Progress Report 1 for basic nomenclature) 


hinge moment modified to effect of axial compres. 
and shear force 
bending moment along a beam strong axis by! 


ree 


“distance | between neutral axis and compression flange 
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Pie: coefficient of dilatation for stress increment in x-direction 


= 


= edge r moment per unit ee to proce unit rotation of edge 
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Plasticity Related to to Design “a 
he This | paper is‘the fourth in a series of reports son 
d Plastic Design emanating from a Joint Committee of 
_ the Welding Research Council and the American 
Society of Civil Engineers. Reference should be 
q to Progress Report No. 1 (which appeared on 
July EM Journal) for such items as the Foreword, = 
"Synopsis, Definitions, and Nomenclature. For infor- ‘ 
mation the complete Table of Contents appears at the 
beginning of Progress Report No. 3 immediately pro- 
herein this paper. Articles 6.3 and 6.4 are included 
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tions the members must have sufficient lateral support so that lateral- a . 


“a torsional buckling will not reduce the resisting moment at these sections. _ 


- When a member buckles laterally it undergoes lateral bending and twisting. | 7 
Elastic regions of the beam resist this buckling through lateral bending stiff- : 
ness ly and torsional resistance. The latter consists of two parts:* ioe 2 
a) ‘Resistance associated with St. s 


A member with an adequate combination of these properties will have effective 
“resistance to lateral-torsional buckling. (6. 29) 


oe ‘The types of lateral support available in actual structures vary greatly. 
For most structures, however, the of also 


points of lateral support. 


Although 1 many on leteral- -torsional buckling of elastic membesi 
have been carried out since Timoshenko (1906, 1913), Wagner (1936), Kappus 
_ (1938) and Goodier (1941), there are few solutions available in the inelastic 
“range. ~ Bleich(6.19) suggests that this inelastic effect can be allowed for by 
_ determining the most critically stressed section and using a tangent modulus 
_ reduction factor (T = Et/E = Gt/G) in the elastic solutions. This approach 
forms a simple and expedient means of approximating the inelastic effects. q 
_ It is generally held that this practice is conservative. However, recent inves- 
; 7" made by Campus and Massonnet(6.32) indicated that in some cases it 
may be on the unsafe side. Clark and Hill (6. 33) refer to these and other ap- 
_ proximations in connection with application to design. In Ref. 6.34 it is shown 
that test results for laterally unsupported wide-flange shapes agree very well 
= the concept that the ratio between the elastic critical buckling stress and 
a actual critical buckling ¢ stress i is the ‘same in the inelastic range and beams 


Neal(6 -35) studied this problem by considering, in detail, the effect of 
(idealized stress - strain curve) o on the lateral f flexural rigidity and 


modulus load, assuming no unloading of the cross section. 
ee Horne(6.37) extended these results for pure bending to the case of I- seams 
His conclusions, however, were modified by other practical considerations. | 


Eventually, from this work, limiting critical lengths of I-beams were pre- = 


=. as a basis for design by Baker, Horne, and Heyman. (6.38) The funda- 
mental concept behind this work is a tacit reliance on the post-buckling 


strength ofa member to provide the necessary moment resistance over a 


ven 


See, example, net. 6. 29 for discussion of st.  Venant’ torsion and 
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the reduced modulus solution of an axially loaded column. | 
ointed this out and deduced an equivalent formula on the basis of atangent 

re 


it inves- 
cases it 
er ap- 


. The lowest point of bifurca 
tion of equilibrium is obtained by considering no strain reversal. At a 


The solution of the various cases on a systematic basis was obtained by 
White{6.31) He also presented a design procedure for determining the spac- 
ing of lateral supports such that the members would not buckle laterally be- fy 
fore the required plastic hinges had formed. _ Although a ‘simplifying approxi- 
mation has been made for practical use, the tavhes still presents three main 3 ; 


1) Under conditions of the moment diagram has a small 
gradient, the procedure gives results which are too 
compared with test results. ah 


ing influence of adjacent segments on the segment being braced. 


3) One of the necessary parameters is the amount of plastic hinge rotation 

=” required to obtain a complete re-distribution of moments corresponding _ 
_ to the ultimate load. Unless an assumption is made of a single given = 
amount of rotation, this would have to be computed analytically oad - 


to torsional deflection takes place before buckling. This approach leads be 


would therefore complicate the design 


_ Assuming as in the solution to local that the 


of yielding of structural steel is non-homogeneous, -30) ‘ the actual material — 
will be in either the elastic or the strain-hardened state. _ Since the material - 

in each of these regions behaves in a linear homogeneous manner charac- — 
terized by its appropriate moduli, it is reasonable to assume that the rea 
can be reduced, ina somewhat simplified manner, to the consideration of a = 
member with an abrupt change in stiffness.(6.27) In the case of a = 


strain-hardened member the resistance to buckling depends almost entirely 

on warping torsion instead of on a combination of St. Venant’s torsion and 
lateral bending strength as in the fully elastic case. For a member with i 
tial strain-hardening, the effects of St. Venant’s torsion, moment gradient _ 
and end restraint are then considered as modifications of the basic solution. i 


_ The basic assumptions for the theoretical ; -analysis(6- 31) ay are as follows: 
| 1) The material in a member deformed beyond the elastic limit either is 


- at a strain less than the yield strain or has reached the strain- — kt 


a) The distribution of the elastic and strain-hardening zones depends on the _ 
magnitude and location of the over- -all inelastic deformation. Strain-— 
hardening is assumed to start at the end of the member at the point of ; 
maximum moment and to spread from that 


3) It is assumed that the member is initially sti straight and that. no lateral or 


this point of buckling, therefore, the linear incremental stress- eee. 


relations can be used, the two different zones being characterized by 


the incremental elastic and strain-hardening moduli, reapectively.° — oe 


*Another solution would be to for the ultimate 
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5) The cross section is ‘plelaeneend to an open peees symmetrical, thin- 


The « — s section is idealized so that the whole section will be con- 
_ ‘sidered as either elastic or strain-hardened under flexure. The condi- 
7 tion of having a non- homogeneous material over the cross section is » 

difficult to handle analytically, | since then the differential equations of 
Pa equilibrium would have variable coefficients. As shown in Ref. 6.31 | 
this: seman of an idealized cross section gives conservative results. 


When the loading is ‘restricted to the c case of loading in one principal plane 
at the boundary only, with zero axial force and zero applied forces along the 
_ beam, the basic equations can be re-written in the following form: miler “ 


Yat” 


u 


bending rigidity of the beam axis 

= bending rigidity of the beam about weak axis a 


- = an le of twist about shear center 
displacement of shear weak 


axis along the beam, 0 <z 
_ The problem 6 solved i in ‘Ref. 6.31 was a singe ~Span assis with load applied 


¥ ~ sumed to be able to resist twisting. Part of the beam was considered elastic 

_ and part strain-hardened as shown in Fig. 6. 16. For the strain- hardened ~ 
ce portion (length « aL) in the above equations were | modified by introducing the : 

strain- hardening moduli in lieu of the elastic values. The condi- 


= 
17) 
— 
Pasi 
Le 


of Strain ‘Hardening 


Dar 


Pook 


= 
6.16. Moment-Curvature Diagram for a Par-— 
9 Yielded Beam Under Moment Gradient 


St. Venant’ 
tioned above were solved bya numerical procedure based on finite differences, 


using the appropriate boundary conditions. The computations for determining | 
the eigenvalue for various conditions of restraint, loading, and extent of , 
strain- ~hardening were carried out on a digital computer. 
series of results were obtained as follows: 


1) Member at each end 


end z simple at end z = 


4) simply at « = 0, fixed at end 
The basic critical buckling for simply beams 


mor 
= moment is from expression A=L. 


the eigenvalue, e equals 11.28 and Z is the plastic modulus of the i section. ae . 
‘For a WF beam it has shown that vi hl y/Z #1. 60 r , therefore 


Ler/ty = = 18. This length is then 1 corrected ee the ‘effects | of moment gradient — rs 


and extent of yielding « over r the uabraced length of the » member. = 


~ 
> 
results) 
| 
a — 
applied 
d fixed: 
as- ia 
elastic 
1g the 


be 


= correction for moment gradient 


factor for partial yielding (aL) 


Vs = factor for St. Venant’ s 


Experimen orrelatic a! ion 
_— specimens ecdititins ali -end conditions by boxing in the end sections 
of the members were tested under constant moment by White. ‘The “fixed-end’ 
condition was one in 1 which the beam could rotate freely in the loading plane 4 


= was prevented with respect to the longitudinal a axis of the beam. _ These 7 
_ conditions, however, were not realized in the actual test as evidenced by _ 
movement of the loading of the end supports. 


- Simple Supports” 
Supports 


— 


g ffect of of the the critic 
— 
— 


2 
sections = 


ced- ent’ 


“Fig. 61 18. Effect Saint t Venant’s s Torsion on the 
nad as: Critical Buckling Length for Beams Simply 8 Sup- 


the values of al obtained from the test results were a oer an 
values predicted by the theoretical analysis for the fixed-end i 
ae Further tests were conducted by R. Sarubbi(6.39) in order to find the effect _ L 
of moment ratio P on the critical buckling length of the simply ~~ ios 
beam. Four specimens were tested with moment ratios of 0.40, 0. -71, 0. 90 -_ 
and 0.98 respectively. Fig. 5. 6. 19 shows a typical curve e of moment v vs. lateral 
rotation for P = 0.71, together with the relation of the moment and rotation =— 


about the strong axis. . The plastic moment My is computed using the average 

value of static yield level | Oys for the flanges and the value of the plastic 
modulus derived on the basis of measured dimensions. _ The non-dimensional . — 
inelastic angle change H/L§, is computed from this curve for three cases: 
the point of buckling, the point of maximum load, and the last recorded load. 

The ratio H/ 1p, is used to express rotation in which H is the total inelastic 7 

angle change between the two load points, L is the distance between the points 

of lateral supports, and Bp = = ‘M,/ EI is the elastic part of the curvature at the — 
plastic 1 moment . The observed 1 values of H/ fy ) are shown on Fig. 6. 20 at the fod 
corresponding ‘moment ratios and slenderness ratios, ry, ¢ of the ‘sections 
tested. These tests that the theoretical conserva- 


000 = lc 
+ 
ig 
Lou buckling length L,,/ry is giver 


| 


TERAL “ROTATION (Rad 
Test Under Moment Gradient == 
tah as mentioned above and these correction factors can be obtained by atrial 
: and error procedure using tables and figures given in Ref. 6.31. For instance, 
q the correction factor vq for the effect of partial strain-hardening of the mem- 
ber depends upon the moment ratio P and the yielded length aL. The value 
me of a depends also on the moment ratio P and the angle of rotation of the “a 4 
beam at hinge points. it may be too tedious to compute the critical length | 
- bettiees two adjoining | points using many figures and tables. Therefore, the 
approximations are for each correction factor. 


4 Then vq can be expressed in terms | of p with a parameter » of rotation 


angle. _ A reasonable value of H/LBp = 3 is chosen on the basis of Pas 


-. following values for WF beams: 
4 


ad web depth. 


<i 
— 
“M, 
| 
ne 
@ 1S Obtained irom the moment ratio and angie Oi rotation H 
4 on the 
4 


Fig. 6. Critical Ratio as a we 
Moment Ratio for Lateral Buckling of Beams 


restraint) are simplified, and the final relationship ese the bee: length : 
and moment ratio can be expressed as a straight line, 


rlty need not be less than 30, since the as is known to be overly 
a The correction factor v. for end restraint depends | upon the a 
of the adjacent segments as shown in Fig. 6.21 where kg is the larger 
stiffness and ks the smaller stiffness of adjacent segments fhe 


cy — (6.39) The stiffness factor k under thrust is approximately é. 


eL is the length of the adjacent segment and is 


__ by Eq. 6.19 for the adjacent segment if it is plastic-elastic or, if the a: 

adjacent segment remains elastic then Ler is given by 


— (1- Aap) - 1100 


| 


— | 
‘ 
— 
of — 
lm 
| — 
— 
= — 


ar 


= 


~ ( for fixed-support ) 


21. Restraint Factor for Continuous Beams rg cau 


= moment ratio in len 
last term gives the of ry to plastic 


for Pao > 0.9, and if <0. 9 vanishes. “res 


btained by entering Fig. 6.21 with kg 
a... ue for Vy is obtained by yen ering ig. wi Ss ky and : 


I euiity is the actual — of the segment being investi- 


Design Approximation for Equal EI/ 
Fora continuous beam with and “equal unbraced ‘segment 


Then ‘Vy is fo found from Fig. 6.21 as before. 


Since the value of ky L/3EI is less than 1. 0, as ‘can aid 


of Fig. 6.21, the variation of Vy ve nearly a = 


AY 


_ 


— 


i 
— 
2.0 
— 
—— 
where Pap = moment ratio in the adjacent seg sig 
int 
: 
— 
| 
— 
— 


(0.9 +key) #2 (oy 


Recognizing that the effect of end restraint, extent of yielding, | St. 
“torsion, moment gradient, required hinge rotation, and axial force are inter- » : 
related and that the exact value of one of these parameters: cannot be deter- 
‘mined without knowledge of the others, a practical won approach would be to 

. assign a representative value to the restraint factor vy i in lieu of the con- - 
sideration of the separate effects. Fig. 6.21 it may be seen that an 
average e value o of i": y is. is 1. 25 for the continuous beam with constant EI and equal, : 


unbraced segment lengths. An equation developed on this basis is as follows: 


ax M 

cr= = (60 - - 40 not less than 

This equation shown in 6. 20 2 as a dashed line. M/M 


v = 1. when 1. 0, » .28. Such a variation is reasonable be- 
cause the correction factor v, for the mpantt ratio M/M increases with ii 
Thus, an approximate 1 value has been established for the effect of all of ~ E 

correction factors required by the theoretical analysis. In the event the de~ 
signer would like to use unbraced lengths greater than those indicated by this pool 
expression he may use the method based on the theoretical analysis by taking oe 
into account each parameter individually and performing the necessary itera- 

tion if he feels the slightly increased allowable unbraced vn les which might 7 
result would be justified by the additional design time aa Oe 


‘The spacing of ene in n the region of plastic hinges ‘aages‘where 
rotation would be required may be determined the 


625 


= end moment ratio 
‘To somewhat ‘greater unbraced lengths, the procedures 
The spacing ng of bracing in the non- -yielded portions of the © 
structure may be methods in * 


”- 
+4 — 
P 
4 
— 
| 
ati 
— 
— 
aw 


The design is a . trial ame error u 

7 braced length is determined by the consideration of lateral-torsional buckling 
: only. However, in practical design the unbraced length may be fixed by other 
structural considerations and then must be examined for lateral-torsional a 


_ buckling. The general procedure for investigating a given unbraced length is: 
‘4 * 1) From the moment diagram f find moment ratios for this segment and its 

(2) Find | Ler/ry fo for each ch segment from Eq. 6. 19 or Eq. 6 1 "aa | 


3) ‘Compute of the adjacent ‘segments from Eq. 6.20. moet 
4) Compute the ratio and the value 
5) Find v,, from Fig. 6. 21. 
+6) The actual length of | of the Segment being investigated sh¢ should be than 


 Dlustrative Examples mals (PF 
n 


the un- 
; wan length that must be investigated in a design will be the one with the A 
largest moment ratio at a hinge point where rotation would be required. Or-— 
dinarily this means that the bracing spacing would be checked in the region | 
of formation of the first plastic hinge. * ‘Or course, , if all segments meet the 
requirement that would be stipulated for a first p plastic hinge location, then 
no further consideration of hinge sequence would needed. 
_ Example I: Continuous Beam as shown in Fig. 6. 22. not notary 
a The example will first be worked without any regard to sequence of hinge — 
_ formation. The problem will then be solved sagen: on the basis that the first 
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Example I of Lateral 


of is discussed in Chapter | 
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BC or CD. The moment ratios are =p = -0. 5, p = — 


2) All spans are in the plastic-elastic state, therefore the critical net 
each span considered as a simple By, = is Gives 

LABer _ LpEcr ae 


rive 


ewe 


which 
ich 


a 

means s that ‘there i is no restraint | on span BC due to. ‘span CD. From 


~BCe: Cc r 


the with the added information that the first 


Since the resulting negative, Kcp = 0. For isons a 
fixed-end (Fig. | 6. 21) v, 1.43. Then 


AGS 
(Lapler = 1-43 x 78.75" = 112.5" 6" (OK) 
Thus, the chosen length is satisfactory. In fact, this was evident in step 2 es 


where it was that (LaBer = 78.75"" with v, = 1.0. 


— 


lt the soproximation of Eq 6.25 had been § 
40 (-0.50) 125-100" (OK) 


a 


Kling 
4 
4 
; 
| — 


Il: (Fig. 6. corner of rectangular frame 


The critical segment is BC. wit. J 


(1 - 0.6% 67) 


= 
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Note: By Eq. 6.25, Cader ace 
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the discussions it has assumed that the structures ‘under 
| consideration were subjected to all of the applied loads acting nner 
It has | been further assumed that the appiied loads increase in their -magni-— 

tude until the limit load is reached. During this proc process none of the applied ~e 
loads are permitted to change their directions and the ratio of loaad magni- ia 

tudes | to each other remains fixed. ‘Such loading is termed proportional load- — : 

a In practice the separate loads on a structure may not satisfy the condition | 
of proportionality. Instead, the separate loads may change in magnitude inde- : 
pendently from each other. These non- proportional loads are called “variable ‘ : 
repeated loads” in the literature. Although these loads act generally in 
random manner, it is convenient to assume a | cyclic pattern in as. the 


ee Fatigue is failure of the material by fracture as a result of repeated load _ 
ing on the structure. In particular, failure will occur at a relatively low _ uy 
lumber of cycles when the repeated loading is such that yielding of the ma- 
terial occurs : alternately in tension and compression at a given cross i. 
This phenomenon is known as “alternating plasticity.” During each cycle of - 
load wit take ‘will lew lead to frac- 


aM 4 | EM 
3) k 
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“ae 
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1959 


b) Deflection Stability: 4" 


: 


This t type of failure is an increase in deflection 
e 


ach cycle of loading, the increments of deflection being in the same direc- | 


tion. _ This mode of failure is termed “incremental collapse.” ‘The problem 7 | 

is to determine the limits on the loads for which these increments — 
= a few cycles of load application and the deflection “stabilizes.” When 
structure reaches this state of stabilized deflection it is said to have an 

: “shaken down” and the corresponding set of loads is referred to as the f 

“stabilizing load” or the “shakedown load.” The st structure re henceforth Lied 


) made by 
47,6. 


~The load P is assumed at first as being applied in a Adeniedicn't direction _ 
the resulting moment- curvature relationship at the | built- -in section is shown 
; from (0) to (a) in Fig. 6.24. If, instead, ‘the load were applied in the opposite - 
direction, the corresponding M-@ curve would be that shown in Fig. 6. 24 
id from (0) to (c). If at point (a) the load P is gradually released and finally 7 

; B ~- in the opposite direction, the M- g relationship is linear for a range 
of moment designated as The magnitude of AMy is l¢ less than or at most 


Ds ia 


= ~ equal to 2My, the precise se value being « a function of the residual stress and = 
Pad the Bauschinger effect. At point (b), yielding starts in the opposite direction. 
Finally a point (d) corresponding to - -Pmax: is reached. To complete the 
rn ae cycle, as loads are released and on reversed, the resulting behavior would 
be as shown by the dashed line d-e-f-a. 
i Failure due to alternating plasticity will x not occur when ranges of moment | 
values exist for which a section behaves elastically regardless of its previous 
loading history. As a first aeperenanien the Bauschinger effect may be 
_ ignored and this range of — ( ris 4 in Fig. 6.24) mes be takenas Tv 
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6. 24. ‘Moment- Curvature Relationship at the Built- in 


of a Cantilever Beam Under Alternating 
where Mj denotes the elastic moment values at any section “i” ” being investi= 
gated. Procedures for calculating the limit for alternating plasticity in 
: of indeterminate structures may be found in Ref. 6.53. atl Me 


ie As noted earlier, under repeated application of a certain esas te loads be 
an increment of plastic deformation in the same sense may occur during each OE 
cycle, The maximum load for which these increments cease after a afew 
cycles is called the stabilizing (or shakedown) load. In Fig. 6.25 is ; plotted ~~ 
_ diagrammatically the the number of loading } cycles versus the deflection under the 
_ load at the end of each cycle. When P is equal to or less than a certain criti- 
cal value Pg, a set of residual moments will be set up in the structure after 
a a few cycles during which the deflection ey ar a limit value. . All er : 


— flection does not stabilize but continues to grow aaa each cycle of load appli a 
: “cation. _ A description of this phenomenon may be foundinthe several 


references th hi a 
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‘Fig. 6. +25. ‘Diagrammatic Representation o of Deflection 


It to mathematically the maximum load Ps for 
4 which the deflection of the structure will finally stabilize. The condition = 
e g sum n of the residual moment My and the elastic moment M ‘produced by | the 
‘ Xe loads does not ¢ exceed the full plastic moment value @ Mp. In general t terms: 
Equation 6.28 may be applied to shown in Fig. 6. 26a, a 
— beam of two equal spans, the supports of which can take upwards and down- 


; q wards reactions. | / Suppose the load P is first applied at B. An elastic analysis > 


- would produce the moment diagram shown in Fig. 6.26b with Mp = 13/64 PL. 
As a second phase of the loading, two equal loads applied at B and D would i 


te give the elastic moment diagram shown in ' Fig. 6.26c. Due to any inelastic ~ 


deformation that may occur under load, the only possible shape of the aikaed co 
residual moment diagram is that shown in Fig. 6.26d. | 
2 ‘ae Eq. 6.28, and defining the residual moment diagram positive as 


in Fig. 6.26d, ‘the to be satisfied are as follows: 


At section B: 13/6 64 (P L) + + | - 

section Cc - -12/64 + 


4 a q 

6.30) | 
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Moment 
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for Moment diagram 


- > for acting at B 


| and (Elastic 
diagram: ofter 
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te 
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However, the residual n moments can have only a — variation across the 


span with the maximum value at C (Fig. 6.26d) sothat 


= 17200 (mM 


pressions stor the stabilizing load P. we obtained from Eqs. 6.29, 6. 30 and 


M 
= 96/19 My/L 06 Mp 


‘The results of this calculation are shown in Fig. 6. 27. The elastic mo- 


ment diagram due to Ps acting at Bis shown in re 6. aa. _The moment ~ 


Figs. 6.27a and 6. 27b, one obtains the final moment Ps act-— 
at B 27d) and at B and D (Fig. 6. 27e). is evident 


‘Thus, from Eqs. 6.32 and 6. 34 


Equation 6.35 indicates that the > stabilizing load is theoretically about 16% 
lower than the ultimate load P, in this example. 
As a check on n alternating the following condition must be 


ang In the previous example at section B ads 


-9/38 Mp | 
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27 -Span Beam of Fig. 6. 26 Under 
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2 \ojaabenttyy, tee increase in deflection for a load P <Pg will cea 
hardening sets in. A discussion of the of 
either in Ref. 48 or Ret. 6.51. ny 


Experimental studies have been performed using continuous beams and — 
frames. Kloppel(6- 54) investigated two-span continuous beams 
with simpie supports. Such structures were also investigated later by _ 
Massonnet(6.50) and Gozum.(6-51) Massonnet tested two-span continuous 
— with loads applied at the center of each span, while Gozum_ investi- J 
gated the same structure with loads off center to simulate | the worst possible 
Neal and Symonds (6. 45) have deflection stability, testing small 
scale rectangular portal frames with symmetrical vertical load and horizon- ; 
_ tal load. Test results indicated that the observed stabilizing loads were oe: 
_ 10% higher than those given by theory. - McCarthy and Popov have experi~ 5 
> 4 - mented on a series of frames under variable repeated loading, (6. 52) investi- 
"gating both alternating plasticity and deflection stability phenomena. Test di 
igi results show that the theoretical analysis (based on no strain- -hardening) © ] 
Gg _ gives conservative predictions and that strain-hardening may increase the © 
All results indicated that the experimental values were higher than those | 
= given by theory. Tables 6.2, 6.3, 6.4, and 6.5, contain a summary of each of — 


previous important observations. In these tables is indicated the loading 


proportio used ina given cyclic load test. The terms a and denote the 


Relationship to Design 


; i problem of variable repeated loading may be disregarded for building frames 
a for the usual conditions Of static loading. The probability of failure” 


dead 

particular is the fact that the ratio of live load to dead 
load must be very large before the load-carrying capacity is reduced because 

* - of load repetitions. In nearly all of the tests described in Section 4, extreme i 


& _ examples were chosen in which all of the > load was « considered to be live load. 
f & is unusual to find such extreme load variations in building structures. x The | 
live load is seldom more than two-thirds total load and uoually it is of 

4 


order of one-third of the total. 
‘Tt must be remembered that the e safety sina does not provide for possi- 
4 mis “ble overloads alone. It also accounts for such additional factors as variation | 
‘material properties, dimensions, workmanship, fabrication, methods of 
q me analysis, etc. Therefore, "variation in live load alone could not ee be 4 
Ss to account for the full value of the factor of safety. = a ; 
Neai(6. 48) has emphasized the fact that failure due to increase of 
is a gradual process so that ample warning of danger is available. This im- 


that a lower load factor 1 is for Ps than provided for 
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Continuous 50) 


Loads applied at mid- -span- 


2 Bent ‘about strong, oxis) 


LOADING 


= Theoretical. ‘load by simple 


the of the most recent rolled have 
| shown that the observed stabilizing load was always greater thanthe 
theoretically predicted value. . Since the theoretical values of Ps are seldom 
_ more than 20% below Py, the practicality of this problem loses much of its © 
Significance . Although account could be taken of variable repeated loading by 
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Deflection stability need not be in the 


design of statically-loaded building frames. 
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oe Laboratory, Department of Civil Engineering, Lehigh University. 
This group has been engaged through the Institute of Research in an investiga- 
tion entitled “Welded Continuous Frames and Their Components. e: The results 
of that investigation were the basis of preliminary reports which were later | 
revised and supplemented by the Joint WRC-ASCE Committee. The — ¢ 


> 


members at Fritz Laboratory included Beedle, G. C. Driscoll, 
_ T. V. Galambos, R. L. Ketter, T. Kusuda, G. C. . Lee, T. Lee, L. W. Lu, 
_ A. Ostapenko and B. Thurlimann. Wm. J. Eney is Director of F Fritz Labora- 
and Head of the Civil Engineering Department. “4 
_ The organizations that supported the research projects out of which | the a | 
_ preliminary reports were prepared were the following: _ American Institute of 
‘Steel Construction, American Iron and Steel Institute, Column Research Coun- 
cil, The Department of the Navy (Office of Naval Research, of 
and pees the Research Council. b 
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PLASTIC DESIGN II 


= Bending moment about strong axis at z= 0 


= Range of moment in which M-@ curve is linear © a 


= Proportion of a given length in strain- hardened 


rf. = Inelastic angle change between points of lateral anal 


Correction factor | for  Venant’ 


Correction factor for extent of yielding 
= Correction factor for end fixity oy 
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Closure by G. — and Thiirlimann 


a 


—Professor Bijlaard is 
“connect See out that the authors’ formula “gon not reduce to the 7 q 
elastic buckling stress when the elastic ‘moduli are substituted. Formula = 
was derived for the strain-hardening range of mild steel as shown in the thesis a 


on which the authors’ reference 3 was based. In the strain-hardening range a 
Professor Bijlaard’s formula gives slightly larger values for the critical 
stress than formula (13). Because of the small value of the tangent modulus, — 
house, the differences in in critical strain values would be considerable. — : — 
Furthermore, Professor Bijlaard’s discussion corrects the statement made 
by the authors that his theory is of the deformation type. Although his theory 
as it is used today in aircraft design is a deformation theory, his complete _ oe : 
theory as published in 1938 includes both deformation and incremental theory. 
cig Professor Bijlaard has been able to reduce all subsequent contributions a 


inelastic plate buckling such as the papers by Dyushin, Stowell, Handelman-— 


_ Concerning the case of the authors’ paper it may be well to quote = 


_ and the one by the authors he discussed to his original theory. aS 


-Bijlaard’s own statement presented at the Annual Meeting of the Column 


Research Council, May 1957 (Ref. (a), p. 55): 


thes writer’s theory (Prof. Bijlaard’s theory) of plastic buckling of x si 


_ plates is in excellent agreement with test results with aluminum plates — 
_ (references 12-15). It is also in good agreement with tests on structural 
7 q steel for small plastic deformations as occur below the yield stress. = 
3 - However, in the horizontal | part of the stress strain diagram at the yield 
bY a stress where the material flows in layers and thus behaves discontinu- 
- ously, in contrast with the continuous behavior assumed in the theory, 
it obviously does not apply* and is too. conservative, as was pointed out 
for columns in reference 16(b). _ This was taken into account by using _ Fe 
requirements b/t< 8.5 and 34 from reference a(c) for hinged flanges 


. and plates where large strains at the yield stress are necessary.” a 
| statement obviously implies that Prof. Bijlaard himself did not ‘asides 


vi that his theory could be extended to the case treated by the authors prior to 

- Instead of starting from a general theory of plasticity the authors determin- xe 
ed expressions for the plate buckling stresses based on orthotropic behavior Be 


of the material. Agreement should exist between these expressions and ‘anol x 


| a. Proc. Paper 1581, April, 1958, by G. Haaijer and B. Thiirlimann. 


1. Research Engr., Appl. Research Lab., U. S. Steei Corp., Monroeville, | a 
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on plastic reduction 
_ advantage of giving general expressions for the buckling strength of ortho- : 


a plates was to clearly indicate the influence of the various stiffness 
factors that are involved. This approach enabled the authors to show en 
source of diffenences in the buckling stresses predicted by various theories | 
_ The authors sincerely — Professor haet s careful review and 
5 
(a) P. Bijlaard, “Buelling of Columns, with Special Attention to 
a Eccentricity, Plasticity and Local Buckling,” Proceedings of the 
ot a Seventh Technical Session, Column Research Council, (May 1958. 
Ss Ref. _(b) P. P. Bijlaard, “Exact Computation for the Buckling of the Webs | 
po Sh Weds} Truss Members Used in Bridge Building, ” De Ingenieur in Ned. Indie 
= No. 2, pp. 1. 17- “1. 23, (1939), 7 


Ref. (c) Ss. Beedle, B. Thirlimann, and } Ketter, “Plastic Design in 


‘Structural Steel,” Lecture Notes, Lehigh University, Summer Course, 
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Closure and Errata by Frank L. DiMaggio 


= FRANK L. DIMAGGIO,! ASCE.—The author v wishes to thank Miss 
Mathews for her pertinent comments. 
the plastic hinges is not known a. priori. . The method presented is valid re- ‘a 
gardless of where the first hinge forms, but the specific results presented in 
the paper obtain only for the case when it forms at the knee. It certainly — a 
would be advantageous to consider the effect of the axial forces on the yield <7 
| criterion in some cases, but the method presented is not amenable to this 
For those cases where many elasto- -plastic stages occur, the ‘author agrees — = 
cat the recent method proposed by H. H. Bleich is far superior. Miss — 
_] states that the author’s procedure is applicable when the load duration is — 
| compared with the period of the first elasto- -plastic mode. . Since the elasto- ig 


plastic response is governed primarily by the rigid body mode when such a 
mode exists, and since this mode has an infinite period, it would be more = 
_ | correct to say that the load duration should be short compared to the time 
to the end of the (first) elasto- ~plastic stage. 
5 -In eq. | (10) change to P; 
- t equation number er for eq. (33) 


P. 16 - - Ine eq. change sin a in first term to: sinl 
interchange expressions a # 


“Ine eq. a. (89) 
an | In eq. (91) change } 3 
a This is not at a typographical | error, but a mistake : in the original - work. te 
ame all numerical results in the elasto-plastic stage incorrect. "y 


18 first paragraph under title First Elastic > Stage ‘change te 


22-In Fig. 7, deflections after first elastic stage are due kod 
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in eq. (91). if correct equation is a maximum deflection | of 665 


“this “exact” This should be expected because the load 345600 


23 - -In ‘Fig. Zo is incorrect errata re. 


- 24 - On the half page starting with Tf, eqs. . (90)” and ending with “ at the top 
of the column,” the numerical results quoted and r — to are incorrect ~ 


: aa Under “Response to Load with Fo = = 57600,” Fi Figs. ,10 and 12 are incorrect 


due to error in eq. (91). Corrected graphs constitute Figs. 1 and 2 of a ; 

errata. The rigid body mode of course still contributes almost all of the — ’ 

: response in the plastic stage. Note that in this case, where the load F_ is only 

a 8 times the static collapse load, a rigid plastic analysis gives results which ‘ 

_are too low. This will always be true in problems of this type because ina _ 
rigid analysis, the coordinate Po has no initial velocity. Ih 


- 9 is incorrect (see errata for p. 22 above) : 
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Ts by Gideon P. willich 


GIDEON P. R. VON WILLICH, 1M. M. ASCE.— Mr. | Chen’s contribution is isa 
Sisk It is noteworthy that the minimum values obtained for the buckling pressure — , 
‘parameter Per by the two methods of minimizing the potential expression do a ihe 
= differ widely. As explained i in | the original paper, the writer feels thet _ 


__ As the writer’s method of the yields v values of p 


strietions on n the mode | of deformation of the shell. 
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a Paper 1897, 1959, by P. R. von Willich. 

, South Africa. 


| 
Ag 
: 
ably be of interest to investigate in detail the mathematical implication of &g | 
. ¥ 
4 4 ‘ 
— 


am. 


— 
— 
— 
— 
— 


"Discussion | by George 


— 


"GEORGE F. ASCE. —The authors have a 


A _ The first of the various test series with which eompuatinen is made is the 
one carried out at Cornell University by the writer and two of his colleagues. 
This information has long since been | published and, therefore, reference _ ie 
should be made not to toa preliminary, eee report (ref. 9) but to a: 
Mason, R. E., Fisher, ( G. P. » and Winter, George: Eecentrically Losted, 
4 Hinged Steel Columns. Proc. ASCE (J. Engr. Mech. Div.) EM 4, Paper - 
1792, October 1958, with discussion in Paper 2099, July 1959. “ogee ba: 

_ It worth noting that in the latter paper comparison is made of the test Fe- 
‘sults: with the SO- ~called R. Cc. interaction equation,” tl the authors’ Eq. 9. 


represented on their Fig. 12. As noted in the quoted paper, the general inter- ae: 


action formula Eq. 9 has been found in satisfactory agreement with tests on a an ‘4 ‘ sd: 


"great variety of cases not covered by the authors’ analysis. Consequently, 
while only “semi-rational,” Eq. 9 has a wide range of validity, having been 7 a 
_ checked by various } investigators against tests on both aluminum and steel 5 i, 
members, in major and minor axis bending, without and with torsional defor- 
| mation, and for a variety of ratios of end moments. 
| Objection is sometimes raised to the presumed complexity” of the inter- a 
~ action formula of Eq. 9 for practical design nuse. In discussing Eq. 9 the pee 
| authors have | shown that the form of this equation is the same as that of their we 
"proposed Eq. 5. Since the American Institute for Steel Construction has seen ; 
‘fit to incorporate Eq. 5 in their Rules for Plastic Design and Fabrication, the a 
complexity of Eq. 5 and, therefore, of Eq. 9 (if suitably transformed) is appar- 
FE not excessive for practical design purposes. 


Naturally, the writer is pleased that the relation which the authors have | a 


found to hold for concentric columns, Eq. 11, is in substantial agreement with ee 

_ that given in another, earlier paper o of which the writer was a co-author (Ref. et _ i 

It is gratifying to note that as a consequence of research effort 

a considerable variety of places the last t ten years, as I reflected call 


a Proc. Paper 1990, April, 1959, by Theodore V. . Galam bos and 


Prof. & Head, Dept. of Structural 
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columns is at last being brought under control. The authors should be 
on hed valuable ‘building block in this edifice. 4 
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